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THE NIMBUS ENERGY BALANCE 


COMPUTER PROGRAM 
ABSTRACT 

The lengthy and laborious hand calculations now associated with the 
design and analysis of satellite power systems can be greatly reduced 
through the use of the Nimbus Energy Balance Computer Program. By 
combining the known electrical characteristics of the power supply com- 
ponents, this program will simulate the operation of three general types 
of power systems - Nimbus, Series Maximum Power Point Tracker and 
Parallel Maximum Power Point Tracker - as the spacecraft passes 
through a complete orbital cycle. The effect of load changes, battery 
failure or changes in operating characteristics, solar array degradation 
or partial failure, and changes in electrical component characteristics on 
power system energy balance can be easily obtained; energy balance is 
defined as the condition where sufficient power exists to supply the 
spacecraft loads and completely recharge the batteries during a given 
orbit. 

The output of the program is a summary of the operating condition 
of the various system components at user-specified time increments dur- 
ing an orbit: solar array maximum available power, actual array power 


supplied, solar array current, battery current and voltage, the relative 
battery state-of-charge and number of ampere-minutes >f capacity re- 
maining in battery, total regulated bus current, peak load current and 
shunt dissipator current. The average power dissipation in the batteries, 
the battery ampere-minute charge-to-discharge ratio achieved during the 
orbit and the battery depth-of-dischargc are also calculated and printed 
out. 


viii 


THE NIMBUS ENERGY BALANCE 


COMPUTER PROGRAM 

I. Introduction 

The design and analysis of satellite power systems is greatly enhanced by 
the capability to simulate the power system's orbital performance. The laborious 
hand calculations usually associated with this simulation can be significantly 
reduced through the use of the Nimbus Energy Balance Computer Program. The 
determination of the various system component operating parameters is ac- 
complished by combining the known electrical characteristics of the solar 
array, battery, source control devices, load power conditioning devices, charge 
controller, system power losses and spacecraft load profiles. A running tally 
of the various power system operating parameters is provided throughout the 
simulated orbit; these parameters are printed out at equal user-specified time 
increments during the orbit. 

This program accommodates an energy balance analysis of three general 
types of power systems: 

1. The Nimbus B series regulator type system (NB), 

2. A parallel maximum power tracker system (PMPT), and 

3. A series maximum power tracker system (SMPT). 

Figure 1 presents a block diagram of the computer model of the three power 
system configurations. By supplying a particular "System Kay" input data card. 
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PTO, VTO (SMPT) 
VAB INB, PMPT) 
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SMPT = SERIES MAXIMUM POWER TRACKER 










the program user specifies which system configuration he wishes to simulate, 
and the computer "switches" shown in Figure 1 are positioned appropriately. 

All input data and computer instructions are supplied by the program user in a 
"data deck" consisting of 50 NCODE cards defining various power system 
parameters, 30 STINT tables containing solar cell, battery, temperature and 
load profile information, and up to 25 Panel Description Cards defining the 
solar array configuration to be simulated. 

Because of the numerous power system components and functions that must 
be considered when simulating orbital performance, the energy balance program 
is divided into a main routine, called MAIN, and the following subroutines: 
STASH, STINT, DRAIN, AMPS and PRINT. All five subroutines are used when 
simulating either the NB, PMPT or SMPT system configurations and perform 
identical functions for all systems. Briefly, the functions of the subroutines 
are: 

MAIN - Perform energy balance iterations 
STASH - Supply solar cell information 
STINT - Store input data tables 
DRAIN - Determine value of load current 
AMPS - Determine available solar array current 
PRINT - Present output data in proper formats 
A more complete description of each subroutine is contained in Section IE. 
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II. Power System Components 


The following paragraphs describe the characteristics of each of the com- 
ponents comprising the computer model of the three types of power systems 
shown in Figure 1. 

1. Solar Array and Isolation Diode 

The solar array may contain fron one to twenty-five solar cell panels con- 
nected electrically in parallel. Each panel may have its own number of series 
and parallel cells, its own solar incidence angle (which must remain constant 
throughout an orbit) and its own temperature -vs -time profile. All panels use 
the same solar cell as the basic building block, and pass current through the 
isolation diode ADIODE. The value of ADIODE (if applicable) includes the diode 
and slipring voltage drop. Total solar array output current, IA, is the sum of 
the individual panel currents at the solar array operating voltage VA. 

The array current -voltage (I-V) curve values are determined by the com- 
puter multiplying the solar cell I-V points by the number of series cells for 
voltage, and by the number of parallel cells for current on each panel. The 
charged-particle-degraded solar cell I-V curve is supplied by the user as an 
input table of I-V pairs along the curve. Temperature coefficients and voltage 
and current degradation factors are also supplied as an input by the user; sub- 
routine STASH manipulates the tabulated cell I-V curve to account for the vari- 
ous design factors and temperature effects. 
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2. Series MPT Unit and Solar Array Bus 

If the program user has elected to simulate the SMPT system, the switches 
in Figure 1 are placed in the SMPT position by the computer and the Series MPT 
Unit is connected in series with the output of the solar array. The component 
transfers solar array power (VA-ADIODE) x IA to the series tracker unit out- 
put with an efficiency, PTEFF, specified by the user. The power tracker output, 
PTO, is then defined by the relationship PTO = (VA-ADIODE) x IA x PTEFF, in 
wauS. PTEFF retains a constant value during an orbit, and a minimum drop of 
1.0 volt is maintained across the series tracker unit at all times. The series 
tracker unit output voltage is designated VTO. 

When either a NB orPMPT system is to be simulated, the switches in 

\ Figure 1 are automatically positioned accordingly; the series tracker unit is 

shorted out, and the voltage on the solar array bus (VAB) is defined as VAB = 

\ 

\ VA - ADIODE. 

3. Shunt Dissipator 

The shunt dissipator is employed in the NB system only, and is represented 
by the equivalent circuit shown below. 


SOLAR ARRAY BUS 
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Whenever the solar array bus voltage (VAB) exceeds the shunt dissipator thresh- 
old voltage (TVSR), shunt dissipator current (ISD) will exist, as determined by 
the effective shunt dissipator resistance (ERSR). The values of TVSR and ERSR 
are specified by the user on input NCODE cards. 

4. Charge Controller 

The NB and SMPT systems use the same charge controller model, shown in 
the figure below. 


IB (AMPS) 



When the voltage drop across the charge controller (VCRR) exceeds the 
dead zone voltage (DZCRR), battery charge current (IB) will increase linearly 
to the point where the maximum permissible charge current (IBMAX) occurs at 
the charge controller knee voltage (VKCRR); further increases in VCRR will 
maintain a constant IB value. If the value of battery charge voltage reaches the 
maximum permissible value (VBMAX), the computer will reduce the charge 
current to a value such that VBMAX is not exceeded, just as in actual voltage 
limiting circuit operation (tapered charge operation). 
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Values of DZCRR, VKCRR, IBM AX and VBMAX are specified by the user 
on input NCOPE cards. 

When a PMPT system is being simulated, the series dissipative charge 
controller previously described is disconnected and replaced by the parallel 
maximum power tracker unit, which transfers power from the solar array bus 
to the tracker output with a user -specified efficiency PTEFF. The tracker out- 
put power, PTO, is normally used to charge the battery, but can also deliver 
power to the peak load and/or the main load regulator if necessary. 

The parallel tracker unit also contains the current-limiting (IBMAX) and 
voltage -limiting (VBMAX) circuits which control battery charge current just as 
in the other two system configurations. In addition, the parallel tracker unit 
continuously compares the total amp -minutes into the battery with the total 
amp-minutes taken out of the battery during the orbit being simulated. When 
the ratio of these parameters reaches the value specified by the program user 
as the C/D ratio (CTOD), battery charge current is automatically reduced to a 
nominally low (0.6A) value. This simulates the actual operation of an ampere- 
hour counter charge control method. If the program user does not wish to take 
advantage of this ampere-hour charge control program feature, a high value for 
CTOD (e.g. 100.0) should be specified in the NCODES. 

5. Battery 

The computer model assumes that all batteries connected in parallel in 
the power system have the same electrical characteristics and can therefore 
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be lumped together into one equivalent battery having the combined capacities, 
the sum of the maximum charge currents and the same voltage as the individual 
batteries. 

The voltage of the computer model of the battery depends upon its state of 
charge (SOC) and the value of current (IB) going into or out of the battery. The 
full capacity of the battery BAMMAX (Battery Amp-Minutes MAXimum) is de- 
fined by the user as an input and the state of charge (SOC) is calculated by the 
computer to be the amp-minutes in the battery at any given time (ACCUM), di- 
vided by BAMMAX. Storage cell data is read into the computer from an input 
data table which tabulates cell voltages as a function of SOC and charge or dis- 
charge current, IB. The computer multiplies the cell voltage by the number of 
series -connected cells in the battery (FUDGE) to obtain the battery voltage VB 
at a given SOC and a given value of IB. One battery is assumed in the computer 
system, having the combined capacity of all the storage modules, and charging 
or discharging at the total value of current. An example of how storage cell data 
is tabulated for computer input is presented later in the report. 

Values of BAMMAX and FUDGE are specified by the user as an input to de- 
fine the fully -charged system storage capacity in ampere minutes and the num- 
ber of series connected cells in the battery. (NOTE: the variable NBAT is also 
used in the program to signify the number of series -connected cells in the 
battery). 
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6. Discharge Diode, Unregulated Bus and Peak Load 


In all three system configurations, the battery discharges through an isola- 
tion diode in the battery discharge path between the battery and the unregulated 
bus. The voltage drop across this diode is called VDIODE, arid is user-specified. 
Current exists through this diode only during battery discharge, providing a 
definition for unregulated bus voltage, VU, during satellite nighttime: VU = VB - 
VDIODE. 

Note that VU can be much greater than VB during solar array illumination, 
since the battery discharge diode is then normally reverse-biased. During 
satellite night, the spacecraft load demand (except for peak load) is determined 
by the computer in subroutine drain to be ILT amperes demand at the unregu- 
lated bus voltage (VU), while ILT is defined at the array bus voltage (VAB) dur- 
ing satellite day. 

The peak load is a user-supplied input data table which defines peak load 
power in watts as a function of time in minutes during an orbit. The table is 
prepared in the same format for all three systems, and in fact can define any 
type of unregulated bus power profile from a minimum duration of one minute 
up to an entire orbit duration load. As seen in Figure 1, the peak load in the 
NB and SMPT systems is supplied from the unregulated bus. The battery will 
discharge if the solar array bus cannot supply enough power at VU to satisfy 
ILT and the peak load current (IPKLD). 

When the PMPT system is selected, the peak load is supplied from a peak 
load bus which is isolated from the unregulated bus by the battery discharge 
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diode VDIODE. This permits the solar array to operate at its maximum-power 
voltage at all times with the PMPT system; the peak load current is supplied by 
the parallel tracker unit up to the limit of IBMAX - any additional current re- 
quirement is supplied by battery discharge to the peak load bus. The value 
VDIODE is applied to the voltage drop across the diode between the battery and 
peak load bus in the PMPT system. 

7. Load Power Conditioning Devices 

Figure 1 shows four types of power conditioning devices which derive un- 
regulated DC power from the unregulated bus and supply the spacecraft loads 
with the desired type of voltage. Any or all of the four devices may be used with 
any of the three systems - the user defines which STINT table location (explained 
in Section IV of this report) contains the load vs time profile of a device which is 
to be employed in the system. The characteristics of each load power condition- 
ing device are described below. 

(a) PWM Voltage Regulator. - This device is a down-converting, switching 
regulator that basically transfers power at a relatively constant per- 
centage efficiency. The iiiput current can be less than the output cur- 
rent in this relatively efficient device. PWM regulator losses (transfer 
efficiency losses) are included along with other system losses in a sepa- 
rate stored table, as a function of regulated load power. The load 
profile for this device is tabulated as spacecraft load watts, at the 
regulated output voltage, versus time; calculation of losses is auto- 
matically made. 
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(b) DC -DC Converter. - This device supplies regulated DC power at an 


output voltage which can be higher than the input (unregulated) voltage. 
The load profile is prepared as regulated output power versus time; 
losses are calculated as a power transfer inefficiency with a user- 
supplied value of converter efficiency EFFCNV. The computer does 
not need the value of output voltage for its calculations; the user must 
define the appropriate constant value of EFFCNV, in percent, for the 
particular device he is simulating. 

(c) Inverter. - This device supplies an A-C output voltage to the load pro- 
file which is tabulated as required watts versus time. A constant per- 
centage power transfer efficiency EFFINV is supplied by the user as 
an input. Values of output voltage, frequency or power factor are not 
required by the computer; the user must supply an appropriate power 
transfer efficiency EFFINV, in percent, for the particular inverter he 
is simulating. 

(d) Series Dissipative Voltage Regulator, - This device supplies current 
to the spacecraft load at a constant regulated DC output voltage which 
is less than the input voltage, similar to the FWM regulator. However, 
in the dissipative device, the input current is assumed to be equal to the 
output current. The product of this current times the voltage drop 
across the regulator is the power lost in supplying the loads with this 
device. The load profile for the series dissipative regulator is 
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tabulated as amperes demand (at regulated voltage) versus time. Values 
of shunt power loss in the series dissipative regulator can be accounted 
for in the system power loss table, explained oelow. 

8. System Power Losses and Fixed Losses 

The system lors watts (SL) shown in the block diagram of Figure 1 are 
values of watts, stored in a table in the computer, that represent the measured 
PWM voltage regulator losses, telemetry and standby circuitry losses, solar 
array bus -to -unregulated bus diode losses and regulated power required by the 
eight charge controllers in the Nimbus B flight power subsystem. These are 
c^lectively called system power losses and are strongly dependent on both the 
unregulated bus voltage and the total spacecraft load demand on the power sub- 
system. Figure 2 presents the measured system power losses as a function of 
regulated bus output power for conditions of satellite nighttime (battery dis- 
charge, low unregulated voltage), solar array illumination (middle-range unregu- 
lated voltage) and shunt dissipator ON (highest unregulated voltage) for a typical 
Nimbus power subsystem. 

When a PWM regulator is specified, the total value of system loss is ob- 
tained from the table which contains file Figure 2 data. If a PWM regulator is 
not used in a system, the losses are calculated as described in 7b, c, and d, and 
in addition the system loss (watts) at G watts PWM load is obtained to account 
for such losses as are caused by telemetry, standby circuitry and other power 
losses associated with the system. 
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Figure 2. Typical Nimbus System Power Loss 
Versus Regulated Bus Load Power 

Figure 1 also shows a computer "switch" which can be positioned to place 
a "power loss, night" (PLN), or "power loss, day" (PLD) load on the unregulated 
bus. The user can specify this constant value of power loss in addition to all 
other loads and losses by supplying the value of watts for daytime and/or night- 
time system operation on an input data card. An example of the use of PLD is a 
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15-watt fixed loss in the series tracker unit of the SMPT system. This fixed 
loss is in addition to the loss caused by the power transfer efficiency of the 
series tracker unit, PTEFF. 

III. Program Description 

This section presents a functional description of the subroutines in the 
energy balance program (MAIN, STASH, STINT, DRAIN, AMPS and PRINT). A 
basic funcaonal block diagram of the computer program is shown in Figure 3, 
which summarizes the important features of each subroutine. 

A. MAIN - Computer Program Control 

The purpose of MAIN is to load input data, initialize power system param- 
eters, select the proper set of energy balance calculations for a particular sys- 
tem, perform a clock function during the orbit, call on the five subroutines for 
data as required and maintain and update values of the system parameters 
throughout the orbit. MAIN employs iterative processes to determine the vari 
ous system voltages and solve for the various branch currents in the power 
system such that the Ohm and Kirchoff criteria are satisfied to within au arbi- 
trarily small error. In addition to the various system voltages and currents, 
MAIN keeps track of battery relative state of charge (SOC), depth of discharge 
(DOD), accumulated ampere-minutes (ACCOM), net power dissipation in the 
battery on an orbit -average basis, ampere-minute C/D ratio achieved during 
the orbit, solar array maximum power and power at the actual operating point, 
and solar array temperature. MAIN also ensures that the maximum values of 
battery charge current and battery voltage are not exceeded. 
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STASH 

DEGRADE SOLAR CELL !-V CURVE FOR 
CURRENT, VOLTAGE AND RESISTANCE 
EFFORT 

STORE CELL l-V CURVES OVER 
TEMPERATURE RANGE 
PROVIDE CELL I OR V FOR A GIVEN V 
OR I AND A TEMPERATURE 


STINT 

STORE ALL INPUT DATA TABLES 

INTERPOLATE IN TABLES AND 
PROVIDE DATA AS REQUESTED 


PERFORM 

SYSTEM 

CALCULATIONS 

CALL 

SUBROUTINES 

MAINTAIN 
TALLY OF 
PARAMETER 
VALUES 


DRAIN 

DETERMINE LOAD AND SHUNT DISS. 
CURRENTS AND SYSTEM LOSSES 


AMPS 

DETERMINE TOTAL SOLAR ARRAY 
CURRENT FROM CELL, PANEL, 
SUN ANGLE, VOLTAGE AND 
TEMPERATURE DATA 


PRINT 

PRINT REQUIRED INPUT DATA AND 
CALCULATED SYSTEM DATA IN 
PROPER OUTPUT FORMATS 


Figure 3. Subroutine Functions in Energy Balance Program 
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B. Subroutine STASH - Solar cell degradation and temperature effects 


The purpose of subroutine STASH is to account for all factors, with the 
exception of charged particle irradiation degradation, which effect an individual 
solar cell's I-V curve. Illumination intensity change, coverglass transmission 
loss, current and voltage measurement uncertainty, ultraviolet effects on cover- 
glass adhesive transmission and standard solar cell accuracy uncertainty are 
examples of factors which are accommodated by translating the solar cell I-V 
characteristic parallel to its voltage and current axes. STASH changes the 
cell I-V curve shape to account for external series resistance and operating 
temperature effects. 

The input to this subroutine is an irradiation -degraded I-V curve*, specifi- 
cation of percentage degradation factors, and the solar cell temperature co- 
efficients. 

After applying the degradation factors to the input solar cell I-V curve, 
subroutine stash expands the degraded cell curve into a family of 15 I-V curves 
spanning the user-specified temperature range of interest for the ensuing 
energy -balance analysis. The 15 I-V curves are held in memory, supplying 
current or voltage data when interrogated during program execution. 

A detailed technical discussion of the techniques employed in the solar 
cell I-V curve manipulation by subroutine STASH is contained in Reference 2. 


*A solar cel I I-V curve, accounting for charged particle irradiation effects, can be obtained by 
employing techniques described in Reference 1. 
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C. Subroutine STINT - Data Storage and Retrival 


Subroutine STINT stores in tables, and then supplies on demand, the values 
of variables which are functions of one, two or three arguments. The first use 
of STINT loads the tabular data into the computer. Subsequent STINT calls 
will ask for a linear interpolation to be performed which corresponds to the 
supplied argument values. 

One and two argument functions may be stored in single tables - Tables 1 
and 6 are examples of one and two argument functions, respectively. Three 
argument functions must be stored in consecutive two argument tables; the 
number of two argument tables must be equal to the number of argumentg values. 

Subroutine STINT is called by means of a statement with the following 
format: 

CALL STINT (ARG 1, ARG 2, ARG 3, FCT, KEY, NGRIPE, MINTBL, 

MAXTBL) 

The first three dummy variables are the three arguments; zeroes must be 
substituted for unused arguments. FCT is the variable value to be calculated, 
and KEY tells the computer what kind of operation is to be performed (data loading 
or data interpolation). A value of -1 indicates the table loading mode, while 
a + 1 calls for a linear interpolation to be performed. NGRIPE is an errc^ flag, 
advising the user of improperly prepared input data, MINTBL and MAXTBL are 
the number of the STINT tables used; both tables are the same if the variable 
is a function of one or two arguments. MINTBL is set equal to the table number 
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which contains the lowest value of argument values and MAXTBL is set to the 

O 

table which contains the highest argument values if the variable is a function 

O 

of three arguments. 

D. Subroutine DRAIN - Total Load Current Calculation 
Subroutine DRAIN determines the total load current demand at the system 
operating voltage for each time increment during the orbital cycle. The sub- 
routine obtains the value of load power (or current in the case of a series dis- 
sipative regulator load) from the profile table for each of the power condition- 
ing devices specified by the user. By adding the system losses and resulting 
inefficiencies to either PLN or PLD, the subroutine computes the total load 
current. For example, during satellite night: 


ILT = 


PWM + SL 
VU 


PCONV 

EFFCONV 

VU 


PINV 

EFFINY 

VU 


PLN 

VU 


+ ISER 


where 


ILT 

VU 

PWM 

SL 

PCONV 

PINV 

PLN 

ISER 


= total load current demand @ VU 
= system unregulated bus voltage 
= PWM regulator load demand, watts 
= system losses, watts 
= converter load power, watts 
= inverter load power, watts 
= fixed nighttime power loss 
= series dissipative regulator current demand 
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EFFINV = power transfer efficiency of inverter 

EFFCONV= power transfer efficiency of converter 

In addition, DRAIN compares the operating voltage with the shunt dissipator 
threshold voltage (TVSR) and calculates the shunt dissipator current (ISD) if 
TVSR has been exceeded. 

E. Subroutine AMPS - Solar Array Current Determination 

Subroutine AMPS determines the total solar array current available at the 
solar array output voltage (VA), accounting for the series -parallel arrange- 
ments, sun angles, panel temperatures and blocking diodes associated with up 
to 25 solar-cell panels. 

F. Subroutine PRINT - Output Data Presentation 

Subroutine PRINT receives the user-specified input data and various calcu- 
lated system parameters from MAIN and writes the oulput tape, preparing the 
data in the proper formats and column headings. 

IV. Program Usage 

The following describes the mechanics of using the program in "non- 
programmer" language. A complete Fortran IV program listing including the 
MAIN routine and subroutines can be found in the Appendix A; Appendix B illus- 
trates a typical data deck set-up. 

The assembly of the complete program as it is submitted to the computer is 
shown in Figure 4. This assembly basically consists of two parts - a program 
deck and a data deck. The program deck, which can be used in either Fortran 
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Figure 4. Assembly of Complete Energy Balance Program and Input Data 
IV or binary form, is always used and is placed first in the assembly. It con- 
tains the MAIN routine and the five subroutines used in the program (DRAIN, 
PRINT, STASH, STINT, AMPS) and does not require any card change from run 
to run to perform its function. 

The data deck contains all the numerical information the program requires 
for computation and defines the user-selected options for each run. Consequently, 
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the data deck must be prepared specifically for each run, or series of chained 
runs, to be made. Cards and tables in the data deck must be positioned in the 
order shown in the program assembly in Figure 4. The input data deck descrip- 
tion and format are presented below in the proper assembly sequence. 

Date Card 

Col. 1-2 Number of month 
Col. 3-4 Number of day 
Col. 5-6 Number of year 

STINT Table Title Card 
Col. 1 Blank 

Col. 2-72 Any alpha-numeric information, such as "DATA TABLES 
FOR NIMBUS D POWER SYSTEM ENERGY BALANCE" 

STINT Tables 

The STINT tables are stacked one behind the other in the data deck in as- 
cending numerical order. The tables listed below are present in the data deck 
in the order shown, for a typical NIMBUS energy balance analysis. 

Table No, Contents 

1 2.0 amp series regulator load 

2 System power loss data 

3 ETA vs. Bat. Temp. Nimbus -B 

4 Relative solar cell current vs. incidence angle 

5 11.4 CRL I-V Curve Unglassed 28 Dec. C AMO 
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Table No, 


Contents 


6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 
17 


3 Mos I-V Curve, Flux is 7.9 Exp 13 
6 Mos I-V Curve, Flux is 1.58 Exp 14 

1 Yr I-V Curve, PHI = 3.16 Exp 14, T = 28 Deg. 

2 Yr I-V Curve, PHI = 6.32 Exp 14, T = 28 Deg, 
400 Watt PKLD Table 




NB SA Temp vs. Time Profile, 612 NM 
25 Dec. C, BOL 
25 Dec. C, 1 Yr. Life 
25 Dec. C, 2 Yr. life 
35 Deg. C, BOL 
35 Deg. C, 1 Yr. Life 
35 Deg. C., 2 Yr. Iife^, 


> Storage Cell Data 


Table No. 


Contents 


18 

19 

20 
21 
22 

23 

24 

25 


PWM Reg 
PWM Reg 
PWM Reg 
PWM Reg 
PWM Reg 
PWM Reg 
PWM Reg 
PWM Reg 


Load 150W 
Load 16 0W 
Load 170W 
Load 18 0W 
Load 190W 
Load 200W 
Load 210W 
Load 22 0W 


No XMTR 
No XMTR 
No XMTR 
No XMTR 
No XMTR 
No XMTR 
No XMTR 
No XMTR 
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Table No. 

Contents 

26 

PWM Reg Load 230W No XMTR 

27 

PWM Reg Load 240W No XMTR 

28 

PWM Reg Load 250W No XMTR 

29 

PWM Reg Load 260W No XMTR 

30 

PWM Reg Load 270W No XMTR 


The maximum number of STINT tables that the program can presently ac- 
commodate is thirty. It is not necessary to fill all 30 table locations in STINT 
if the data is not needed. 

The first card of each STINT table is a header card, which must be pre- 
pared in the following format: 

Cols. 1-8: Any alpha-numeric characters can be used for a date. 

Cols. 9-12: Table number. Cannot be zero. Fixed point and right -justified. 

Cols. 13-14: Number of argument^ values. Cannot be zero. Fixed point 
and right -justified. 

Cols. 15-16: Number of argument values. Cannot be zero, is 1 for a 
function of one argument. Fixed point and right-justified. 

Cols. 17-19: Not used. 

Cols. 10-70: Any alpha-numeric characters desired. Usually used for table title. 

Cols. 71-72: 00 

After the header card, each card in the table uses 10 fields of 7 columns 
each for the argument values and the function values. The first card contains 
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the first nine argument ^ values in fields 2 through 10. In the following cards, 
field 1 contains an argument value, and fields 2 through 10 contain correspond- 
ing function values. After all the argument values have been spanned, the whole 
series of argument cards followed by argument cards can repeat until all the 

I Z 

function values are used. If there is an argument value for the table, it goes 

O 

into field 1 of the argument ^ card. Columns 71 and 72 on each card must contain 
a sequence number, starting with 01 for the first card. Figure 5 shows a typical 
STINT table coding sheet for a single argument (solar cell current as a function 
of voltage) STINT table. Figure 6 shows a typical two-argument STINT table 
format. 

After the last STINT table in the data deck, there is a card labled END OF 
STINT TABLES, starting in Col. 21. Cols 9-12 and 71-72 must be left blank on 
this card. 

Run Label Card 

Following the END OF STINT TABLES cr > d is a card containing any desired 
alpha-numeric information in Cols. 2-72, which usually describes the first run 
to be made, such as: RUN NO. 1, NIMBUS D, 1 YR IN ORBIT, 25 DEG C BAT- 
TERY, CONSTANT 150W REG BUS LOAD. 

NCODE 

Following the Run Label Card are the 50 NCODE cards. The card number, 
or NCODE, is right justified against Col. 3. The numerical value of the NCODE 
variable is left-justified against Col. 5 and must have a decimal point. Table 1 


24 



25 


Figure 5. Format for o Typical Single-Argument Stint Table and Header Card 



FORTRAN TABLES FORMAT 
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Figure 6. Format for a Typical Two-Argument Stint Table ond Header Card 




Table 1 

NCODE Names, Numbers, Typical Values and Description 












Table 1 (Continued) 

Numbers, Typical Values and Description (Sheet 2 of 2) 





















shows the NCODE names, the NCODE numbers, the NCODE values and a brief 
description of each NCODE for a sample computer run. Only the NCODE num- 
ber and its numerical value are punched on the NCODE cards; the other data in 
Table I is for information only. All 50 of the NCODES are initially loaded into 
memory, thus a single run or the first of a series of chained runs must contain 
all the NCODES in the data deck. The 50 NCODE cards used in the data deck are 
described below: 

1. TN is the total orbit nighttime in minutes. 

2. TO is the total orbit period in minutes. 

3. AT is the time increment between energy balance calculations in minutes. 

4. VBMAX is the maximum battery voltage permitted during power system 
operation. 

5. IBM AX is the maximum allowable battery charge current. 

6. DZCRR is the dead zone voltage of the battery charge controller. No 
battery charge current will flow unless the voltage across the charge 
controller exceeds this value. 

7. VKCRR is the battery charger controller knee voltage. The maximum 
allowable charge current will flow into the battery when the voltage 
across the charge controller meets or exceeds this value. 

8. CTOD is the battery ampere-minute charge-to-dis charge ratio (C/D). 

The NCODE is needed only when simulating operation of an ampere- 
hour controlled PMPT System, and reduces the battery charge current 
to 0.6 amperes when the prescribed C/D ratio has been achieved. 
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9. TVSR is the shunt dissipator tum-on voltage (volts). 

10. ERSR is the equivalent resistance of the shunt dissipator (ohms). 

11. PLN is a constant nighttime loss, watts. 

12. PLD is a constant daytime loss, watts. 

13. VDIODE is the battery discharge diode voltage drop (volts). 

14. BAMMAX is the ampere-minute capacity of the fully charged battery. 

15. ETA is the panel normal to sun vector angle. A value of 0.0 must al- 
ways be inputed for this NCODE (accounting for panel angle is ex- 
plained later). 

16. NSLT is the number of the STINT table which includes the power system 
losses. 

17. NBMINT is the STINT table number with the minimum battery temper- 
ature information. 

18. NBMAXT is the STINT table number with the maximum battery tem- 
perature information. 

19. NCELLT is the number of the STINT table which contains the input 
solar cell I-V curve. 

20. NPKLD is the number of the STINT table with the peak load power profile. 

21. NPWM is the number of the STINT table with the PWM regulator load 
power profile. 

22. NINV is the number of the STINT table with the inverter load power 
profile. 
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23. NCNV is the number of the STINT table with the converter load power 
profile. 

24. NSER is the number of the STINT table with the series dissipative 
regulator load current profile. 

25. SYSTEM KEY tells the computei* to simulate a series maximum power 
tracker (-1.0), a Nimbus B system (0.0), or a parallel maximum power 
tracker (1.0). If this card is omitted from the data deck, the Nimbus B 
system will be simulated. The NCODE is called PMPT in the MAIN 
program listing. 

26. NPRINT is the time increment in minutes between calculations at which 
output is printed. 

27. NEND is an "end of runs" key. A value of 0.0 implies that additional 
computer runs are to follow, while a 1.0 signifies that this is the last 
run to be made. 

28. EFFINV is the power transfer efficiency of the inverter, in percent. 

29. EFFCNV is the power transfer efficiency of the coi.. ?rter, in percent. 

30. NBTEMP gives the number of the STINT table with the battery tem- 
perature. A value of 3.0 must be inputed for this NCODE. 

31. NDGRAD must be set to 1.0 in the first run. This causes the machine 
to automatically degrade the solar cell and expand it for temperature 
in Subroutine STASH as specified by the degradation and temperature 
parameters in the NCODES. When chaining additional runs, if the solar 
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cell degradations are not changed, NDGRAD should be set to 0.0 in the 
second run. By setting NDGRAD to 0.0, needless repetitive computa- 
tions in the solar cell subroutine are eliminated. 

32. SIGISC is the solar cell short circuit current temperature coefficient, 
(AMPS/°C). 

33. SIGVOC is the solar cell open-circuit voltage temperature coefficient 
and is punched on the card as a positive number; the program later 
gives it the proper negative sign (volts/°C). 

34. Dll is the first short circuit current degradation factor; it usually refers 
to a standard cell error. (This parameter is given as a percentage of 
remaining current after correction is made). 

35. DI2 is the second short circuit current degradation factor; it usually 
refers to a solar illumination intensity variation, (%). 

36. DI3 is the third short circuit current degradation factor; it usually re- 
fers to an ultraviolet degradation (%). 

37. DI4 is the fourth short circuit current degradation factor; it usually 
refers to a current measurement error (%). 

38. DV1 is the first maximum power point voltage degradation factor; it 
usually refers to an external series wiring loss and is specified as a 
percentage of voltage remaining after correction is made. 

39. DV2 is the second maximuni power point voltage degradation factor; 
it usually refers to a thermal cycling degradation loss (%). 
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40. AVFMO is the maximum power point voltage of the undegraded solar 
cell in volts. 

41. AIPMO is the maximum power point current of the undegraded solar 
cell in amperes. 

42. AVOCO is the open circuit voltage of the undegraded solar cell in volts. 

43. THETA is the open circuit voltage degradation factor, it usually refers 
to a voltage measurement error. This parameter is specified as a 
percentage of voltage remaining after correction is made* 

44. TNOT is the input solar cell reference temperature in degrees centigrade. 

45. NBAT is the number of series storage cells in the battery. 

46. ADIODE is the array blocking diode voltage drop, specified by the 
user. A value of 0.0 is punched if no diode drop is desired. 

47. PTEFF is the power transfer efficiency of the maximum power track- 
ing unit, either series or parallel. This parameter is defined as a 
percentage. 

48. DELTT is the temperature increment between the 15 stored solar cell 
I-V curve in STASH (°C). 

49. ADDT is the temperature increment to be added to TNOT to determine 
the highest STASH temperature (°C). 

50. DENFAC is a curve shape temperature correction factor for use in 
STASH. This factor should be specified as 0.065 for a 1 ohm-cm solar 
cell, and 0.0 for a 10 ohm-cm cell. 
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Array Signal Card 


Immediately following NCODE 50 card must be a card containing 999 in 
columns 1-3. This card tells the computer that solar array information is to 
follow. 

NPANEL Card 

Following the Array Signal Card is the NPANEL card, which contains the 
number (NPANEL) of solar panels in the array (maximum number of panels is 
25). This number must appear right-justified in columns 1-3; no decimal point 
is required. 

Panel Description Cards 

Following the NPANEL card is a panel description card for each solar 
panel in the array, up to a maximum of 25 panels. The number of these cards 
must agree with the value of NPANEL. Each card contains four fields of ten 
columns each, in floating point format (requires decimal point). 


Columns 

Variable 

Typical Vaiue 

1-10 

No. of Series Solar Cells per String 

94.0 

11-20 

No. of Parallel Strings per Panel 

36.0 

21-30 

Solar Incidence Angle (degrees) 

0.0 

31-40 

Panel Temperature vs Time Table 

11.0 


Location in STINT 



Following the Panel Description Cards is a blank card. This tells the 
computer to stop reading in data and to start computing. If it is desired to 
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chain an additional run, a new Run Label Card and only those NCODES and Panel 
Description Cards that contain changed or new information should be placed after 
the blank card. In addition, NCODE 27 must be set to 0.0 fo * all except the last 
run, when it must have a value of 1.0. As many runs as are desired can be chained 
in this manner, ensuring that sach new run starts with a Run Label Card and ends 
with a blank card. Refer again to Figure 4 for the proper sequence of card posi- 
tions for chained runs. 

As indicated in Figure 4, computer control cards are required in front of 
the Program Deck, in front of the Date Card, and behind the last blank card at 
the end of the Data Deck. The particular control cards needed vary from com- 
puter to computer, and are sometimes different for identical machines at two 
separate facilities. 

A complete listing of the input data deck, including typical control cards, 

20 example STINT tables, and example run label cards and NCODES for two 
chained program checkout runs can be found in the attached Appendix B. 

V. Program Output 

The information that the computer equipment prints out after an energy 
balance run consists of the following items: 

1. STINT Table Summary 

The STINT table number, the date on the STINT table header card and the 
title of the tabulated data as it appears on the header card are listed for each 
table stored in STINT. 
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2. Input Data Page 

• Run number and date 

• Run comments (as specified on input card) 

• Listing of NCODE numbers, names and values 

• Solar Array description: panel number, number of series solar cells, 
number of parallel solar cell strings, solar illumination incidence angle and 
number of STINT table which contains the array temperature-time profile. 

3. Subroutine STASH Printout 

• Values of temperature for which the degraded solar cell I-V curve has 
been prepared are listed in a row across the page. 

• Values of degraded solar-cell maximum-power voltage and current, 
open-circuit voltage and short-circuit current appear in columns under each 
temperature. 

• Values of every other calculated current and voltage pair comprising 
the I-V curve and stored in the computer memory are listed in columns under 
each temperature. 

4. Power Subsystem Data 

The names of the calculated power system parameters are listed in a row 
across the top of the page: orbit time at which the calculation was made (TIME), 
number of ampere -minutes in the battery (ACCUM), relative state of charge of 
the battery (STATE), output voltage of the series tracker unit in the SMPT 
system (VTO), unregulated bus voltage during satellite nite or solar array bus 
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voltage during satellite day (VU(Night), VAB(DAY)), solar array current (IA), 
solar array power at the operating point (PA), solar array maximum power 
(PM AX), battery voltage (VB), battery current (IB), load current including sys- 
tem losses (IL), shunt dissipator current (ISD), peak load current (IPKLD), and 
solar array temperature (TEMP). 

5. Battery Data Summary 

The depth of discharge, in percent of capacity at beginning of run, the 
ampere-minute C/D ratio actually achieved during the run, the charge energy 
into the battery during the run in watt-minutes, the discharge energy out of the 
battery during the run in watt-minutes, and the orbit-average power dissipated 
in the battery during the run, in watts, are presented. 

If several runs are made at one time, the output data for each run is printed 
out in the same format as for the first run, except that only changed values of 
NCODES are printed out on the input data page, and subroutine STASH printout 
will not appear if the same solar cell I-V curve and degradation factors as in 
the previous run are used. 

Figures 7, 8 and 9 show a typical energy balance computer run output. 
Figure 7 shows the NCODE and array description card listing. A portion of a 
typical STASH printout is shown in Figure 8, while Figure 9 shows the output 
of the power subsystem data. 


37 


RUN NO. 1 ON THIS DATC OF 12-13-6° 


TABLE CTMMrKTS. 

NIMBUS TABLES 

RUN COMMENTS. 

RUN NO 1 

NE V OR CHANGSC PARA*tTERi> 


t TNINIQHT TIME) 75.9900 

? TOIOPBIT T I MF )•«... AO. 00''i' 

3 DELTA 1 .9909 

« VBNAX (VOLTS). s«.0«On 

5 IBMAX (AMPS)... E.^O^O 

6 DZCBR (VOLTS) 0.««"0 

7 VKCRR (VOLTS) 1.A99" 

a c/o ratio... l.noo 

9 T VSR ( VCLTS ) • . • SS.oonn 

10 ERSP (OHMS) O.m ’0 

11 0.0 

12 PLO O.C 

13 VD10OE( VOLTS)....... 0.5009 

1* HAMMAX (A-m) 21 AC .0000 

15 PTA-(OEC) . 0.0 

16 NSL T 2 

IT ISBM INT ..... 12 

10 NSMAXT 1? 

19 NCELLT.. S 

20 NPKLC.. 10 

21 NPWM TO 

22 N IN V. ................................ . 19 

23 NCNV. *0 

24 1 

25 SYSTEM KEY (-1 ,0*5MPT .0.0»NB .1 .0*«»MPT) 0.0 

26 NPRTNT 2 

25 t'FINV as. 0000 

29 EFFCNV 90.0000 

30 NBTEMP 3 

31 NOGRAD 1 

32 S1GISC (A/OFG 0.00014 

33 SIGVOC ( V/DEG C>... 0.00’70 

34 Dll (PERCENT) 9«.0«*90 

35 O I ?(PERCFM I .......... 100.9090 

36 O I 3 (PERCENT ) 100.9990 

37 DI4(PEPCFKT) lO0.0«09 

30 DV1 (PEPCENT ) 99.0000 

39 DV7(PFRCENT ) . ICO. 0009 

40 AVRM0( VOLTS) . . 0. 47000 

4| A |PMO( AMPS) .. 0.17090 

42 AVOCOt VOLTS >•• n.5«99n 

A3 THETA (PtPCEM 100.9990 

44 TNO) T (DEG 75. "900 

45 NO. OF CELLS/PATT. 24 

46 ADInDF( VOLT) , . 1.0990 

47 PTFFF.... 97.9999 

40 PEL TT ( DEG C) 10. t'900 

49 AOOT (DEG Cl...... 37.0090 

50 OENFAC ( 1 9»HM-r<4*0.0 6L.10 UMN-CM*0 .0 ) . . 9.O6SO 

PANEL t 

CELLS * 192. c STRINGS * 132.00 INCDN ANGLF * 0.0 TFMP T A9L c = 1 


Figure 7. Program Output-NCODE Listing and Solar Array Description 
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EDIT OF TrxpEBUTUBE AND DEGRADATION CDPRFF T EO SOLAR CFLL I -V CURVFS FDR FNERGr RALANfF PGM. |N AUTOMATIC IKiriF 
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Figure 8. p rogram Output— Subroutine STASH 









tint accum state via wy ( n ioht ) i* 

(STS UNLT> VABIOAV) 


PI PV4* 


»L 


146. IPKL" TCMP 


0 • D 2160.00 
2.0 21*6.02 
6.0 2131.00 
6.0 2119.4’ 
w.-o Tm.ii 

10.0 2092.1* 

12.0 2070.46 

14.0 206*. 40 

16.0 2030.69 

10.0 2036.72 

70.0 2022.60 

27.0 2000.36 

24.0 19°4.’6 

26.0 IC 42.31 

20.0 2004.9* 

30.0 2017.34 

32.0 2030. 

34.0 20*3.31 

36.0 2036.0? 

36.0 207 C.. 3* 

40.0 2004. Cf 

42.0 2097.07 

64.0 2112.0*3 

46.0 7126.40 

40.0 2140.04 

30.0 2153.30 

32.0 2170.01 

94.0 7104.00 
3STO 2197.7? 

36.0 2209.91 

90.0 2219.66 

62.0 22?7 . 77 
6 ».0 2239.10 

66.0 2»4*,P« 
• 0.0 2*30.00 

70.0 2294.*? 

72.0 2294.91 

74.0 7273.?? 

76.0 2200. 7C 

70.0 2204.00 

60.0 2 293. *7 


1.000 C.C 

C *994 0*0 

t.«p7 O.C 

C.901 w.O 

6.573 0.0 

0.900 0.0 

0.062 O.C 

0.964 C.O 

0.945 0.3 

0.9*3 O.C 

C.c* O.C 

C.Q30 C.O 

C.O?’ C.C 

C.°? 3 O.C 

C.9?« 0.0 

C. 07 * O.C 

C.o* a O.C 

0.9*6 0.0 

r. 03 ? C.C 

r.5*« C.C, 

0.964 O.C 

•07* C.C 

C.97« C.C 

<*.«n* c.3 

0. 901 O.C 

C.9Q0 C.O 

1 0O« O.C 

t.O*? 0.3 

1.0’* O.C 

!.o?3 C.C 

1.007 C.c 

t.C’l C. 1 * 

1.0** O.c 

1. c»0 O.i 

O.C 

1.0*6 C.C 

’.0*5 C.C 

I • "*' ? C.C 

i.">*e c.o 

1.0*« O.C 

1.063 


0.0 

32 . ’*? 
’2. I *■* 
37. ?6 ’ 
?1 . °4’ 
31 . * 
31 . ‘A* 
31 . **» 
31 • ’0° 

n . 1 

30 .96* 

30.796 

30 . **r 

33.9” 
34. « 46 



?*•*’* 
3* . >"* 
35.1** 


35.4’ ’ 
15."" 
35. *4* 


3 7. ’of 
35."’ ^ 
30 



39.0’* 




39 

3? • * ’ r 
38 • " 

3’ . *’*' 



DEPTH 06 Of SCh*RCF(PFBC''»4T ) -<.272 


* . n»6 
. u 1 




, a 0 




C/0 RATIO 1 .702 


CHARGE FNSRC.Vf HATT-Wf NOTES) 


* ««»»r .n** - 


OISCHAPGF EN7RGV ( *ATT -XlNUTtS > 


’OWB AVG PWR CISSIPATEC IN BATTERIES! WATT* ) 6 -.. 7^.4 


Figure 9. Program Output-Power System Data 


VI. Energy Balance Calculations 

At the start of an energy balance computer run, the values of the NCODES 
and other internally -used system parameters are initialized. The data deck is 
then read in: STINT table data is stored in the computer memory, the NCODES 
for the new run are read in, updating the initialized values, and the solar array 
configuration, sun angles and temperature profile are defined from the panel 
description cards. This preliminary effort is done by MAIN, from the beginning 
of the routine through instruction 105 (Refer to Appendix A for listing of computer 
instructions). 


40 


A. Solar Cell Data 


After the initialization described above, subroutine STASH obtains the user- 
supplied solar cell I-V curve from STINT and degrades the curve with the cur- 
rent, voltage and series resistance factors specified in the NCODES, This de- 
graded I-V curve is then corrected for temperature effects and finally 15 I-V 
curves are stored in memory at 15 different temperature values, over the tem- 
perature range (TNOT + ADDT) to (TNOT + ADDT -15 DELTT), as specified by 
le NCODES. STASH will interpolate between the 15 stored I-V curves when 
lied upon by MAIN or subroutine AMPS for particular solar cell data. 

B. Spa Dec raft Nighttime Calculations 

Energy balance calculations for spacecraft nighttime are made in MAIN, 
from instructions 106 through 3005. Since all three system configurations (NB, 
PMPT and SMPT) operate in the same manner during solar array eclipse periods, 
the same set of instructions in the computer is used for all three systems during 
the nighttime duration, TN. 

Initialized parameters for the nighttime calculations are as foUows: 

(a) Orbit time is set equal to zero: T = 0. 

(b) Total ampere minutes removed from battery is set equal to zero: 

TOTOUT = 0. 

(c) The total energy, in watt -minutes, removed from battery is set equal to 
zero: EOUT = 0. 
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(d) The battery relative state-of-charge is set equal to 1.0: SOC = 1.0. 

(e) The total ampere minutes presently in battery is set equal to the initial 
battery capacity, in ampere -minutes: ACCUM = BAMMAX. 

The nighttime portion of the orbit (TN minutes duration) comes first, with the 


following computations being made at each DELTAT time increment. 

(1) Increment T = T + DELTAT. If new T is greater than TN, nighttime 
calculations are completed; if not, continue with step (2). 

(2) Assume a battery discharge current of 1 ma; determine storage cell 
discharge voltage: VBDCH @ SOC and IB = - 0.001A. 

(3) The unregulated bus voltage during nighttime (VA) and battery voltage 
(VB) are calculated: 

VB = VBDCH x FUDGE 


VA = VB - VDIODE 


(4) The total load (ILT) and peak power currents (IPKLD) are calculated 


at VA: 
ILT 


PINV PCONV 
FWM + SL EFFINV EFFCNV 
VA + VA + VA 


PLN 

VA 


+ ISER 


IPaLD 


PKLD 

VA 


(5) The battery table is reentered at IB = - ILT - IPKLD and a new value 


of VBDCH is determined at SOC and IB. New values for VB and VA 


are calculated as in step (3). 

(6) Step (4) is repeated at the new unregulated bus voltage. The sum of 
ILT and IPKLD are compared with IB from step (5). 
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(7) If the difference between the two values is less than 10 ma, go to step 

(8). If the difference is greater than 10 ma, return to step (5) to 
determine a new value for VBDCH at IB = - ILT - IPKLD. 

(8) The ampere -minutes removed from battery during time interval DELTAT 
is computed: 

AMPMIN = IB x DELTAT 

(9) The number of ampere -minutes remaining in battery is determined: 

ACCUM = AC CUM + AMPMIN 

(10) The new battery state -of-charge is calculated: 

_ _ACCUM 
BAMMAX 

(11) The total amper-minutes removed from the battery is computed: 

TOTOT = TOTOUT - AMPMIN 

(12) The energy, in watt-minutes, removed from battery is determined: 

EOUT = EOUT + AMPMIN x VB 

(13) The time is compared with the total nighttime; if T < TN, go to (1). 

If T - TN, the battery depth of discharge is calculated: 

DEPTH = (1.0 - DC) x 100 

Return to step 1; if T and TN have different values, go to step 1. 

When, in step 1, T is greater than TN the spacecraft has emerged into the 
sunlight. The solar array is illuminated and a much more complicated set of 
computations must be made at each time increment. 
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C. Spacecraft Daytime Calculations 


Energy balance calculations for spacecraft daytime (solar array is illuminated 
and orbit time is greater than TN) are made in MAIN, with a different set of com- 
puter instructions used for each of the three system configurations. 

Before the daytime calculations are begun, the value of equivalent charge 
controller resistance is calculated: RCRR = (VKCRR-DZCRR)/IBMAX. In- 
structions 130 through 806 are common to all three systems and are executed 
at each time increment during the daytime before the energy balance computa- 
tions are started; the solar array maximum power PMAX at the orbit time T is 
calculated. The values of total ampere -minutes into the battery (TOTIN) and 
total energy into the battery (EIN) are set to 0.0 before the daytime computation 
begins. 

1. Nimbus B (NB) Daytime Energy Balance Calcuations 

The NB system uses the 900 series of instructions in MAIN for daytime 
energy balance calculations. 

(1) The storage cell open -circuit voltage (VBO) is determined at SOC 
and IB = 0. The battery voltage (VB) is calculated: 

VB = VBO x FUDGE 

(2) The total load (ILT) and peak load (IPKLD) currents are calculated 
at an array bus voltage (VAB) set equal to VB: 


ILT 


PWM + SL 
VAB 


PINV 
EFFINV 
VAB + 


PCONV 

EFFCNV PLD 
VAB + VAB 


+ ISER 
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IPKLD = 


PKLD 

VAB - VDIODE 


(3) Subroutine AMPS is entered at array voltage VA = VAB + ADIODE 
to determine the solar array current (LA) available at the array 
operating voltage, VA. The available battery charge current is 
calculated: IB - IA - ILT - IPKLD. 

(a) If IB = 0, go to step (4) 

(b) If IB < 0, go to step (5) 

(c) If IB > 0,go to step (6). 

(4) The solar array operating power (PA) is calculated: PA = IA x 
(VAB + ADIODE). The battery parameters are updated: 


AMPMIN = IB x DELTAT 


ACCUM = ACCUM + AMPMIN 


SOC 


ACCUM 

BAMMAX 


EOUT = EOUT - AMPMIN XVB 


TOTOUT = TOTOUT - AMPMIN 


Time is incremented (T = T + DELTAT), and program operation 
operation is returned to Step (1) to begin calculations at the next 
time increment in orbit. 

(5) The storage cell discharge voltage, VBDCH, is located at SOC 

and IB and the battery voltage determined: VB = VBDCH x FUDGE. 
The array bus voltage is set equal to VB; ILT and IPKLD are 
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calculated as in step (2). The total array current is determined 
VA = VAB + ADIODE and the actual value of battery discharge 
current is calculated: IB = IA - ILT - IPKLD. Go to step (4). 

(6) The battery charge current is set equal to IBMAX and the storage 
cell voltage (VBM) is determined at IB and SOC. Battery voltage 
is calculated: VB = VBM x FUDGE, and array bus voltage is de- 
termined: VAB = VB + DZCRR + (RCRR x IBMAX). The total 
load and peak load currents are calculated as in (2). Subroutine 
AMPS is entered (at VAB + ADIODE) to determine the solar array 
current. The available battery current is calculated: IB = IA - 
ILT - IPKLD. If the available battery charge current is greater 
than IBMAX, go to (a); if IB < IBMAX, go to (b). 

(a) A voltage increment is defined: DELTAV = 0.1 volt. The ar- 
ray bus voltage VAB is set equal to TVSR + DELTAV and the 
array current at VA = VAB + ADIODE is determined. Values 
for ILT and IPKLD are calculated as in (2). The available 
battery charge current is calculated: IB = IA - ILT - IPKLD. 
(Note that ILT now contains some value of shunt dissipator 
current ISD - defined as (VAB - TVSR)/ERSR. The new value 
of IB is again compared with IBMAX; if IB is still greater 
than IBMAX, set DELTAV = DELTAV/2.0. A new array bus 
voltage is calculated. (Note that VAB is now equal to 
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TVSR+ 0.1 + 0.05). System currents are again evaluated, and a new 
DELTAV is added to VAB in order to further increase ISD and 
reduce IB if IB still exceeds IBMAX. A maximum of 10 iter- 
ations is allowed in this manner to zero-in on the values of 
system currents. Go to (7). 

(b) If IB < IBMAX, the array bus voltage is calculated: VAtj-VBOx 
FUDGE +DZCRR. ILT, IPKLD and IA are determined at VAB, VAB- 
VDIODE, and VAB +ADIODE respectively. The storage cell voltage 
VBMis located in STINT at SOC and IB= IA-ILT - IPKLD; the 
battery voltage is determined as VB = VBM x FUDGE. The 
array current is determined at VA = VB + DZCRR + IB x RCRR) 

+ ADIODE; ILT and IPKLD are calculated at new VAB = VA - 
ADIODE. The new charge current (IBN) is determined; 

IBN = IA - ILT - IPKLD. If | IBN - IB | >10 ma, a new value 
of storage cell voltage VB is located at IB = IBN; this interation 
is repeated until the difference between IBN and (LA - ILT - 
IPKLD) is less than 10 ma. When |IBN — IB j £ 10 ma, go to 
( 7 ). 

(7) The battery voltage is now compared with the maximum al- 
lowable voltage (VBMAX). If VB < VBMAX, go to (8); if VB> 
VBMAX, the computer moves the system operating point out 
along the I-V curve just as in (6(a)). Less input power is 
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obtained at each new voltage point until both the currents and 
voltages are compatible with acceptable small errors and 
maximum limits; then proceed to step (3). 

(8) The battery parameters are updated: 

AMPMIN = IB x DELTAT 
ACCUM = ACCUM + AMPMIN 
EIN = EIN + AMPMIN x VB 
SOC = ACCUM/BAMMAX 
TOTIN = TOTIN + AMPMIN 


The solar array operating power is defined: 

PA - (VAB + ADIODE) x IA 
The final value shunt dissipator current is determined: 


ISD 


VAB - TVSR 
ERSR 


(9) Orbit time in incremented: T = T + DELTAT. If T < TO, go to 
(1); if T > TO, go to step (10). 

(10) The ampere minute charge-to-discharge ratio is calculated: 

RATIO = TOTIN/TOTOUT. Orbital average power dissipated in 
the battery is determined: 


PAVG = 


(EIN + EOUT) 
TO 


The energy balance run for NB is now complete. 
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2. Parallel Maximum Power Tracker (PMPT) Daytime Energy 
Balance Calculations 

The PMPT system uses the 600 series of instructions in MAIN for daytime 
energy balance calculations. Another parameter is also used with this system 
and is computed at each time increment: RATIO = TOTIN/TOTOUT. TVSR is 
set to 1000 volts to avoid unintentional use of a shunt dissipator in this system, 
which can operate at a solar array bus voltage up to 80 volts or greater. Refer 
to Figure 1 for the PMPT system block diagram. 

(1) The unregulated bus voltage is initially defined as the previously- 
calculated solar array maximum-power voltage minus the array diode 
drop, VU = AVMPSA - ADIODE, and the solar array current is set 
equal to the current found earlier at the maximum power point, IA = 

IAM. The total load current ILT is calcu’ated at VU; ILT is compared 
with IA. If (IA - ILT)< 0, go to (2); if (IA - ILT) > 0, go to step (3). 

(2) As the total load current is greater than the array current, battery 
discharge power is needed to support the load. Cell discharge voltage 
VBDCH is looked up by SUNT at SOC and at IB = -1.0 ampere. Un- 
regulated bus voltage is calculated VU - VBDCH x FUDGE - VDIOD5. 

Solar array current IA is found at VU, as well as ILT and IPKLD. Total 
value of battery discharge is calculated: IB = IA - ILT - IPKLD, solar 
array operating power is determined, PA = (VU + ADIODE) x IA, and 
the battery parameters are updated just as in step (4) of the NB day- 
time calculations. Orbit time is incremented and step (1) above is repeated. 
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(3) If (I A - ILT) from (1) is positive, the value of RATIO is compared 
with the user-specified CTOD; if RATIO is greater, the battery has 
reached a full state of charge and a limit value of charge current is 
set: ILIM = 0.6 A. If RATIO is less than CTOD, set ILIM = IBM AX. 

(4) The parallel tracker unit output power is calculated: PTO = VU x 
(IA - ILT) x PTEFF. Battery open-circuit voltage at SOC and 0.0 
amps is calculated and an initial value of charge current is found: 

IB = (PTO/VB) - PKLD/ (VB-VDIODE). Note that this value of IB will 
high since a low VB (open -circuit voltage) /as assumed. If IB < 0, 
go to (5); if IB > 0, go to (6). 

'5) If IB from (4) is negative (this can only occur if a high peal: load 

exists), storage cell discharge voltage VBDCH is located in STINT at 
SOC and IB = - 1.0 A. The battery voltage is defined (V3 - VBDCH x 
FUDGE), and a final value of battery discharge current is determined: 
IB = PTO/VB - PKLD/(VB-VDIODE). The remaining system param- 
eters are defined: VA = VU, PA = (VA + ADIODE) x IA. The battery 
parameters (SOC, AMPMIN, ACCUM, TOTOUT, and EOUT) are up- 
dated just as in step (4) of the NB daytime calculations. The time is 
incremented (T - T + DELTAT) and the calculations in step (1) are 
begun again. 

(6) if IB from (4) is positive, it is compared with the value of ILIM. If IB 
is greater than ILIM, battery voltage VB is found at SOC and at 
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IB = IUM and IB is redefined: IB - PTO/VB - PKLD/(VB-VDIODE). 

If IB is still greater than ILIM, a new operating voltage is specified: 

VU = VU + 0.1. IA is found by AMPS at the new VU, DRAIN supplies the 
new value of ILT and step (3) is repeated. Eventually, the increasing 
value of VU will reduce the available IB below ILIM, as follows: 

Battery voltage is found at the new value of IB, and now another value 
of IB is obtained from VB: IB - PTO/VB - PKLD/(VB-VDIODE). 

This value of IB is now compared with ILIM; if IB is greater, VU ,'s 
again incremented, further reducing IB. After sufficiently increasing 
the array operating voltage, IB will be equal to or just slightly less 
t han ILIM, and the last value of VB is compared with VBMAX. 

(7) If VB is greater than the maximum permissible VBMAX, VU is incre- 
mented to a higher value just as it was when IB was too great. Charge 
current is further reduced at the higher array operating voltage and 
eventually the reduced value ol IB will result in the VBMAX limit not 
being exceeded. At this time the battery parameters (AMPMIN, RATIO, 
EIN, TOTIN, ACCUM and SOC) are updated, array operating power is 
calculated, the values of all system parameters sent to PRINT and 
orbit time is again incremented. 

(8) If orbit time T is less than or equal to TO, step (1) is repeated; if time 
T is greater than orbit duration TO, the actual C/D ratio achieved is 
calculated: RATIO - TOTIN/ TOTOUT, battery power dissipated as 
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heat is calculated: PD = (EIN + EOUT)/TO and the energy balance 
run is complete. 

An additional feature of the PMPT coding instructions is that when a final 
system operating point has been determined from energy balance considerations, 
a check is made to see if the unregulated bus voltage VU is at least one volt 
greater than the battery voltage, during system charge. This check is necessary 
since an operating point is assumed at the maximum -power voltage, which con- 
ceivably could be even less than battery voltage for an unusual combination of 
high array temperature, low battery temperature and a severe solar array volt- 
age degradation. If VU is too low, the computer will increment VU = VU + 0.1 
until a satisfactory operating point is reached, and an error message will be 
printed out a.. the completion of the run. 

3. Series Maximum Power Tracker (SMPT) Daytime Energy Balance 
Calcuations 

The SMPT system uses the 700 series of instructions in MAIN for daytime 
energy balance calculations. As seen in the system block diagram of Figure 1, 
the operation of the SMPT system is identical to that of the NB system after the 
series tracker unit has processed the solar array power. 

(1) The tracker output power is determined at the array maximum power 
point (IA = I AM and VAB = AVMPSA - ADIODE): PTO - IA x VAB x 
PTEFF. Assuming zero current into the battery, the storage cell 
voltage (VBO) is obtained from STINT, and the battery voltage is 
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defined: VB = VBO x FUDGE. The series tracker output voltage (VTO) 
is set equal to VB, and the total load current calculated: 


ILT 


PINY 

PWM + SL EFFINV 
VTO + VTO 


PCONV 

EFFCNV 

VTO 


PLD 

VTO 


+ ISER. 


A first value of battery charge current is defined: 

IB = PTO/VB - PKLD/(VB-VDIODE) - ILT. 

(2) If IB from (1) is negative, battery power is required to support the 
spacecraft load. Storage cell voltage (VBDCH) is determined at SOC 
and IB = - 1.0A and battery voltage is defined: VB = VBDCH x FUDGE. 
DRAIN is entered to determine the new ILT at the lower VB, and a re- 
fined value of IB is calculated as in (1). Battery parameters are up- 
dated, just as in step (4) of the NB daytime calculation, and calcula- 
tion at the next time increment is begun in step (1) after the values of 
system parameters have been sent to PRINT. 

(3) If IB from (I ) is positive, IB is set equal to IBMAX. After storage 
cell voltage (VBM) is determined at IB and SOC, the battery voltage 
is calculated: VB = VBM x FUDGE. An initial value of tracker out- 
put voltage is defined as VTO = VB + VKCRR. A new value of ILT is 
calculated at this value of VTO, as in (1), and a value of IB is determined: 


IB = PTO/VTO - ILT - PKLD/ (VTO - VDIODE). 


(4) If the value of IB from (3) is greater than IBMAX, the system operating 
voltage on the array is incremented: VAB = VAB + 0.1, where 
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originally VAB = AVMPSA - ADIODE. A new value of LA is obtained 
from AMPS, and step (1) is executed again. This procedure is repeated 
until IB is less than or equal to 1BMAX. 

(5) After an IB which is less than or equal to IBMAX has been determined, 
either from (3) or after the last iteration of step (4), a maximum of 10 
iterations are started to find the pi. ope r value of IB. The charge con- 
troller resistance is calculated as RCRR = (VKCRR - DZCRR)/IBMAX 
and the tracker output voltage is set equal to VB + DZCRR, where VB 
is defined as the battery open circuit voltage (VBO x FUDGE). A value 
of ILT is found at VTO as in (1) and a value of IB is calculated: IB = 
PTO/VTO - ILT - PKLD/ (VTO-VDIODE). The storage cell voltage at 
the new value of IB is located in STINT and a new battery voltage cal- 
culated; VTO is redefined as VB + DZCRR + (IB x RCRR). ILT is re- 
calculated at the new VTO, and a new value of IB is defined: 

IBN = PTO/VTO - ILT - PKLD/ (VTO-VDIODE) 

(6) The new value of battery charge current is compared with the older 
value; if | IBN - IBj <10 ma, the error is deemed small enough and 
program operation is shifted to step (7). If the difference is not within 
the allowable 10 ma tolerance. IB is set equal to IBN. The calculations 
of step (5) are repeated, up to a maximum of 10 times, each time re- 
placing IB with IBN and allowing the system operating point to be de- 
termined with little error. 
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(7) After IB from (6) has been determined, the resulting battery voltage 
VB is compared with VBMAX. If VB is greater than the user-specified 
VBMAX, the array bus voltage is incremented VAB= VAB +0.1, as in step 
(4), until the decreased available power at the higher VAB produces a 
low enough value of IB such that VBMAX is not exceeded. 

(8) System operating voltage is checked, just as with the PMPT system, to 
ensure that VAB is greater than VB by at least 1.0 volt. When this 
criterion is satisfied, battery parameters (SOC, AMPMIN, TOTTN, 

EIN, ACCUM) are updated, solar array power is calculated: PA = IA x 
(VAB + ADIODE), and step (1) is repeated for the next time increment. 

If T > TO, RATIO and PD are calculated as with the other two systems, 
and the energy balance run is completed. 

VII. Summary 

A computer program has been described which will simulate the operation 
of three general types of satellite power supply configurations (Nimbus B, Paral- 
lel Maximum Power Tracker, and Series Maximum Power Tracker) as the space- 
craft passes through a complete orbital cycle. The simulation is accomplished 
by combining the known electrical characteristics of the various system com- 
ponents. The system’s operating condition (array power, load powers, battery 
voltage and current, etc.) is determined at each time increment during the orbit 
and printed out. The program will enable the user to quickly and easily ascer- 
tain the effect on the power supply operation of load changes, disconnecting one 
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or more batteries, occur ranee of a partial solar array failure or expected array 
degradation and differing component electrical characteristics. Power dissipation 
in the storage models under any conditions can be determined, as can the battery 
ampere-minute charge-to-Jischarge ratio achieved during the orbit. Excellent 
agreement has been demonstrated between the computer calculated power sys- 
tem operation and actual hardware measurements. 
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APPEltolX A 


-C**** 1NIMBUS POWER SUBSYSTEM ENERGY BALANCE" PRUGRAM 

DIMENSION NUAYE ( 3) 

DI ME NS I O N KS I 2 3) 

DIMENSION XC(23> »XSt 23) , NTtMPN 1 23 > , ANGL It 23 1» ANGLA 
LOGICAL im I TE 


01 0001 0 
01 00040 

frt-tTtrtnrr 

( 2D ) t PTEMPS (25 ) 01 00040 


COMMON NPANEL — “ 

COMMON ATEMP»AVMPfAX IMPURE YE 


READ - < JA»700U)NDATt( 1 > ♦ NUATE < 2) , NUATE t3 ) 

NKUN=1 ^ 


XS<I)=U.U 

XC( I )=o.o- — 


01 00070 
01 ntTTNtT 
01 UOOVO 
Ui 0010 U 
Ul 00110 
01 0 U 12 U 
Ol ooi3o 


ftlMIBilli 


Ol 00130 
01 00140 


ANGL1 ( I 1=0.0 
RS ( 11=0.0 


MUUEL=0v0 
DIl=lUU. 
— tH 2 =-t 3 tn — 
D I 3= 1 00. 
014=100. 


ol 001*0 
— e-t ooi vo 


ill -00210 
Ol" 00220 
o l - tw fr so - 
-01-00240 
oi m >230 


UV2=lOU.O 
Tn=o. 0 


01 002/0 
01 0 O 2 MO 


DEL » aT=0.0 
X I Brt AX = 2U • 0 


Ol 00300- 

01 00310 


VKCKK=1 .40 
CLSK=O.U 
IVSR - lOOO. 
tKSK=0.0 
PL N = U • 0 


VUlUUE=O.U 

BAmMAX=0.0 

Arrti f K+« EDO. 

NttUTU = 1 • U 
NDGKAU=0 
S I 0 I SC = U. 0 


Ol 00330 

Ol 00340 
01 003 3 0 
01 00340 
Ol 00370 

-O l 00 3 * 10 - 

01 003*0 
01 00400 

Ol- -ItO^r lO- 

Ol 00420 
01 00430 
01 0044U 
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SIGVUOO.U 

TnijT=o.o 

» t ) 

NVttMAX=0 
Vt)MAX = 0.0 

X1CGL = 0tU 

fc - 1 ou • u 
PlKfM = 0.0 

P L K O - O. O 

MKLUAO=U 
fcl A = 0 » 0 

WtUA I J « U 

NC H A I N- 0 
NPK IN I =i 

w ) . 

NPL UAU = 0 
tSA I t1*lP-0, U 

— r Oj lfcHP lp 

At) H)0t = 0.0 
Ott_AMG = B*0.0 

YlH 

Tt)YUUT=0.01 
UtG I 0 I ~ 1 • U 

OLNr *00«07 * 0 

Nbt-T=0 

ntp”)n=0 



n)CNV = U 
tvbPK=0 

L - H- {t* V — i -O O. -t* 

tPPCMV= 1 00 • 0 

x iso=o.o 

bY b K P Y - O r- O 

PYfcFP-lOO.O 

CTOl>“ 10.0 

N PKLO - 0 

PL AG=U.U 

fc 10=0.0 

p4H44*O-»-0 * — 

PU=U.O 

C Kt AL) I'AhLfcb CUf..ttlMl LAKU 

itHi H P aD 700 H 

■C LUAtJ SY I inY lAttltb 

CALL b I 1 IM ! ( O. U, O. Ot O. Ot 0.0*-l T HGK I Pt * U t 0 ) 

Il-(WGKIPt) 1UL, lUl, 1 Ok: 


04 -00*-00- 

Ol 00**0- 

01 00 * 70" 

- Ol UMBO 

Ol 00*40 

01 00000 

01 OOMO 

— Ol 00*20 

Ol 00 *- j » U 

Ol 00*5*0 

01 00**0 

01 00 * * U 

Ol 00*5 70 
01 00*HO 

Ol 00 *4 0- 

01 00*00 
Ol 00*10 

01 00*2 0 - 

01 00**0 
01 00*70 

01 OU * H O 

Ol 00*00 
Ol 00700 

04 — 0 0 * 4 ( f 

01 007 ko 
01 00 /BO 

01 007 * 0 ' 

Ol 007*0 
01 007*0 

— 0 4- " 0OY-/ O - 

Ol 00 7 HO 
o4 007*0 

04- OO H OO - 

01 00* 1 0 
01 00*20 

04 — OtA - 40 — 


01 

oon*n 

^ I 

ol 

IM.M i J W 

OOH*() 

01 

OOH/O 

01 

UOHHO 
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102 WRITE IJT},7001) 
WRITE ( JB * 7002 ) 


00R90 

00900 


CALL EXIT 

READ RUIN COMMENTS 


UCHU , O M » f U U J I 

FKIN1 HEADER 

WRITE ( Jb, 7004 INRUN, (NUATE< J) , J»1 , '3I 


00920 
009 TO 


Itllk'Llll 


00950 
009 Ml 


WHITE ( Jb » 70U5 ) 
WKi'fE ( Jb , 7003 ) 


009811 

00990 


1U3 READ < JA,70U7)NCUIJE,FAKAm 

IE (nCUUE-999) 11U3,99,1103 

— l 'rtfs- !> — tn cu t j r ) io 5 , 1U5, r»4 

99 KEAU { JA, 7UU7 iNFANbf 

IE (NFANEL.GT.25) GU “Ttf 4UU 


KEAU < JA,66U0)F1 ,y*,F3,F5 
66UU EOKMAI (10El0.5> 


0101 U 
01020 


01040 


01U6U 
010 10 


XM I )=F2 
anglai i i i 


X 1 I I — ' 

NTEWFS ( I I =F5+.Ul 

6UUU WKlft ( Jb,660l ) I »XC( I ) ,X5( l > , F3 , IV I tMFb ( I t 


15H INCUN AWtitE * E5.1r~ 
14H IEHF lAttLE = I2T 


0 lt)90 

01 TOO 


6602 FU.'PIAT ( 1 HO ) 

GU TU 103 


01 

tJl 


01120 

01130 


01150 

01160 


Oil 80 
01190 


123,24, 25.26, 27, 28 .29, 30, 3 1,32 .33 ,34 .35, 36 ,37, 3 8, 39,40* *rt,42,' til -01210 

— 243,44,45,46,47,48,49,50 TT^t/trOtO! 01220 


WRITE < Jb.HUOl INLUUE, Tr< 
GO TO 103 


01 

01 


t! I 24() 
01250 


WRITE < Jb,8002)NCUUE,Tt.l 
GU TO 103 


012/0 

01280 


WRITE ( Jtt,H003)N(;ilUE,0ELl AT 
GU TlJ 103 
4 VttElAX = FAR Al“l 


01300 

01310 

(1320 
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60 TU 103 
l#> NSLT=FAKAn+.0l 

wa n t? 

GO TO 103 

17 W»r.lNT«FAKAM+.Ol 


01 Ol<M*0 
01 01700 

ill a l 7'n 

V/ 2 ' • i ' y 

01 01720 
01 01730 


GU TO 103 

18 NHWflXTsO AKA-i+.Ol 

WKlTfc ( Jtt, 8018 HCOOb ,NHMAX T 


01 0 1 7 •>{> 
01 01 7 M> 
01 01770 






H *t>6RAtt=RARAr»+.Ol — itl 02 1M) 


■ 80 lit- 103 e 

32 _ Sl'Gf5€ = RARAi' , i — ■ — * i 01 021H0 


GO TU 1U3 01 02200 

33 SIGV()C = RARAM 01 02210 

WKITfc ( Jb,8U33)i\lCUUt, SIGVUC Ul 02220 
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GO TO tt>3 
34 DI1»RARAM 


GO TO 103 
33 UI2«RARam 


- 50 TO 103 

ai. oi i.auiAi* 






j\j ui j-rwnwn 


GO TO 103 

53-014-RARah 


GO T O 103 

A i fci 1 — A 1 1 A ^ _ _________ .. 

jJj UV 1 *rftK«n 


GO TO 103 


—0 2320 


— 01—02360 


GO TO 303 
-42 AVOCO*OAKAM- 


voeoi-» avoco 

GOTO 103 


Oft I TO f 8043 ) NGOUt t T HfcT A 
THETA-THETA/100.0 


01-03530 

ul uCy“v 


T * OT = O ARAH — 

OR ITfc- 4 00 r 8044 INCUOEtTNU! 


01-02560 

01-02570 


43 NBAT=PaRam+O.01 

WRITE < Jtt,8043>NCU0t*NHAf 


50 TO 103 

46 AOlOOt*ftARAm 


-51 02300 
Ol 02600 

l\ 1 A n f 1 


—51 02620 
5102630 


- GO TO 103 

y. ~J niLCC-DA UAm* 


01 

n l 

02630 

/ i ? A 

Hf r ttrr - rfiKAfl 

WRITE ( JB,8047)NC00EtPTEEE 


U I 
01 

02670 








PTEFF=PTtFF/10U.0 
GU TU 1U3 


Ot 02680“ 
Ol 02690 


WRIT t ( JB» 8U48 INCUUt » tJ fc L I I 
GU TU 103 


WHITE ( JB,804y)NCUUt,AUUl 
GU TO 103 


01 02710 

-trt- 02720 


w X V* c_ » .*u 

— 01- 02740 

01 02750 


WRITt I JB , 8050 INCUUt , OttvP AC 
GU TU 103 


til- 02770- 

0 1 02 7 80" 


_ . 

L nUU b C C t r i INI* rULLUWb 

UI u^nuu 

105 I F ( NOGRAU I 50 1 * 5U1 , 5U0 

01"028 IO 


uv=ovi noo.*nv2/ioov 

NOG R A 0=0 


-01-02830 


CAit sTfrsmtnsc," 

lAUUTvOttrfil 


lOtrLTT r 


502 WRITE < JB, 8997t~ 
-- NGRiRt =7777 


rrTivtrf y^tr 


mmimmgmm mihmmmgmmmmmm 

i Kiiriki ici'iviviiiiii nirxiiciki ■^^law^aTimn jAinTiuiu 



— RCRP^t9R-CRR-0ZCKR 1 /X I BMAX 
— +P — ONCttftlN) lOO, 106, 107 


-IP tt)B TEMPI 1109,1110,1109 

L109 T/Att -STiNT tPTA,0.0,0.0,BATfcinP,l,NG-KlPt,NBTtPi8TNBT-fcinPt- 



0 1 03030 
0 1 0 3 0 4 0- 
-0-1-03000 


1110 ntalct=o 

NP«Gtr=50 


— oiro3tr7o 
—01 03080 


CALL PRINT (NPRINT,NPAGt,NTALLY,T,ACCUh,5>UCrOi0,0.0,O.0,0.tr,0.0,0.0l 03100 
10, U. 0,0. 0,0. 0,0. 0,0. 0,0.0) 01 03110 

108 T =T+UELT AT 01 03120 
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6999 CONTIMJt 0 3 031 30- 

IR <VKCRR> I0il3»i0ll2, 10113 01 03390 


30333 

UCLTAV-VKCRR 


01 03360- 


I ftR a 20 


Vl uJl 0X7 
A) A317A 


bli lU iuiin 








iTcR-10 


trr l/3i'Tt7 

10119 

Lufvt iNUfc 


"Ul UJCUtT 


rT trtr j 1 1 ) fc U* U 




RTtmPS < 2 ) =0. 0 


UI ~\JDCCXJ 

309- 

irt T-tN H 0 V 9 1 U V t 
CALL STINT <SOC» 

X i vJ 

-O I 00 3t BA T GwB j VBOCH-; 3 t*»G8 IBB t NBH t W T t NB H A X T ) 

01 03290 


V T 0-0. 0 


-VU0Lrt=VBOCH*R006fc 

IH N GKHM E Ull f 11 H 112 

Hi OGR 1 Rb=l 

— 1 1 3-0*333 ( JB ♦ 7008 >oGR H*t 
GU TO 200 


V B a VBUC H 


N GR I rb< 


CfttL O k fl fft-f XILtXf L T-rRtrNt V^C t i iRt TVSRt lrRirRT^Trt G R l H: ! 

% tiru.-i _■ t „.u -■<• .in ■» . r~ i n l k- i ' l nr r »..n T i.it . 

iNrTftt trri tyV t Nt»r it fr»jTK 9 trrrl n>r ftTTtnV t t rtr 




CALL ' S- TI W T < Tr Oi-GT O» Gf 1 » N GR 1 1‘t: t H l‘ K LU , N OKLO ) 

XltBLLO-RKLO/VA : 

- 137 X1B— X3LT-X1RKLO 

3000 C A LL STINT (SUC T Xi BTB^ T bw H T V B UC Hy l f NGKl^ teTWBH l H T TNBw OXT > 

IF 1NGRIB6-) 30O3 T 3OOl f 3002 

— 3002 NGHlffc ■ 20 - 

GO 4U — H 3 

3001 VbDCH = VBL)CH*hUOGt — — 

VA ■ VHOCH - VO lUOt 

N URI B L ■ r O UJ«U : 


03 03230 

01 032 6 0 

— 01-032 VO- 
— 03 032 BO 
01 032 9 0 

nl n ‘j ra f\/t 

vJT Ujj\7u 

— ot- 03310 

01 03320 

— 01 03330 

01 033 9 0 

0 1 033 6 0 - 

03 03330- 

01 0 3 380 

- -03- 033 90 
-03-03400- 
01 036 1 U 

01 03 44 0 


- Vo ■ VOUCH 

CALL ORAiW (XIL»XILS,BLH t VA,CLbKtTVSK,tKSRtT»HGRlFt»fiSt-Tr 

lN f‘W n f N 1NV T H C lTVt KtitK) t O IN V i h GH C B Vt t H ) 

i 3003 IMXILS + RRLU/Va+XIb-U.OI) 3 U 05 , 3006 , 3006 

3006 XlB=— XILS— PKLU/VA 

G U II J 30U0 

3006 X 1 LT = XI Lb 
TtriRT » 0.0 

T fc nHb " 0»0 

VU = 0.0 
X I A=U.U 
RA = U .0 


— Ol — 0 34 6 0 

r» i / k n y / a 

01 U3 9 7U 

—03 03460 
01 03490 

01 03 6 00 

03 03610 
— 01-03620 

01 03 6 30 

Ol 03640 
01 03660 
01 03660 


i 






I 
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10T4A TtWf> T «PT1?W P S m - - 
6O"f0~tt4 


60 Ttr-t^A 


Cum-gCCtrw-fflmHwtrr 


TOr^rAf rtJT*f t Tt to rAo^ 


01 O3AZ0 
01 "03 630 


— 01 03 6*50 
Ot 173660 


Ol 03680 

01-03 6vo- 


— tit— 03?oti 
01 037 Hr 




^BSBiSBSSwSnSI^BSS! 


1 fOl lit = ~Tt)T I rj*fli*iKm lo 
Knf It) = TOT In / foTUUT- 


v/u I >-/ 1 wu 

- iro-iPtl^Ttrri30,t3ov2oo - 

207 HiKON=ivfU)rfAl 


— or 03310- 



l c.7 

LHLL C A 1 1 

U l u^nnu 

1 BU 

OO 80A"tKt»t=l »NPANtt 








Kt Yfc = U 



tJt tTB^V 

tr 

IP IXCILKrnM 80l»8U2»801 



tTl t)3VO 

0 * 




802 CONTINUE 
KtYE=l 


CA.L STASH (AXIMP, 
1 »UELTT ) 


l i — v ■ v 

00 806 L=1 tNPANtL 

806 PAR=PAR+XS<t> 


PMSAtt= XI AM * < AVMPSA - AOHlOE > 

PMAX ■ XI AM * AVMPSA 


O*** NlMbUS 4i WITH PEAK UJAO REGULATOR STTERTG 

900-VTU-O.0 


- IP < XC < L > I 911*911*912 
912 CALC ST TNT 1 0 * E4 A » 0 *0 1 PTEmPG ft4* TrrNGRO 


— — — &AtL ST4 nT f S0CTU.O-»ttAT-EMP, VbU, 


I WPWM * N I MV »nCnV rNSPft*EEE l N V rtfMCNV *4 N4- 


lNCELLT,AOIO0E*0E-LANGrR5*M0UOtt* VtJCUl *i 


IE <XIb) 9Ul» 902 ,903 

901 CALL STINT < SUCtXItttttATErMPt VBOt 1 tNGRIPEtMOMhvTtNPMAXT I 


— 04— 00 



01 04-290 
Ol 04260- 


V Ab-Vb 01 04280 

CALL OKAIn (XILtXILTtPLO, VAbtCLSHr4VSKfEKS«tTtNGKIPEtiNSLTr— (Hh 04290 


£AtL AMPS- fXI A»VA9fPTEmPST*S?*CtETAf ANGLAt AnGLI fNGRIPbtWAUTOr Ol 04340 

— INC ELL T? AOlOOt t ObL ANG » KS t MUOOtL t VOCO I t AOOT t OEL T T 1 01 04320 


GO TU 904 
902 X I H a 0 • 0 


04 04340 
Ol 043S0 


903 XIri“XI«MAX 

CALt STINT ISOCtXIbtHAfEMPtVbMfitNGRIPEtNbMlNTtNbnAXT I 


VAb=Vb+XIb*KCKR+OZCRR 

CALL ORAIn (XILtXILTtPLUtVAbtCLSKtTVSktEKSRtTtNGKIPEtNSLT, 
INPPM tNl imV f NCNV t NSER » EEE I IV V t tEECN V t T iM ) 


01 04370 
Ol 043 HO 

t . 1 t\l. 'j i t / 1 


Ol 04400 
01 04410 
01 04420 
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CALL AMPS (XrA,VAB,PItMPStX3,XCfETA,ANGrA,ANGLr»NnRIPfc,r«mcrrmr - 0r-r>4V30 
1NCELLT t AUI LJUfc * UfcL ANG t RS t MUUUfcL » VUCtJI t ADHT « OfcLT T) ~ 011)4^(1 



VAb = VAb-tJfcLV 
GO TU 926 
925 VAb=VAB+UbLV 


07 



926 CUiMTlNWfcr 

CALL UK A I N< X 1 1, X I LT t HLU» V AB » CLSK f FVSKf tKSK* 


CALL AMKS<XI«tVAB,KTtnKS,XS»XCf feTA t -AiM6LA,ArrCttri 
1 NC feLL T * AU I OOfc , UtLAwG » KS » «UUUt:L t VUCtri , Afltth t)tt T+ 



CALL ST I NT < SUC » X I B , BA T-fcrrF t VB , 1 , rtGKi Ft , 
VB=VB*PUOGt 


GO TO 924 

fTFKATttm TU FTWO GMONf 

OlSSIKATOK COKKt-itT T-tlttthfS 

906- UtrL f Av=0. 1 

ni'w.ii 



VAB~TVSTOH>tirFAV 

91© X7SO-XH.T-XIC — 



01 04 68 0 

-01 0469© - 


-IF < 1 TtK— 1G1 9l6,916 f 9©4 

916-GAtL ArtFS- t XI A, 9ABrF7£*FS^ XS»*C »fc7A T *f*Gt*T4NGti 


C Att GRAIN <XTL^XHT^FLOtVABfCLSi 

1 k i l~\ . l u ^ i T i 1 i * /* kl Li — i L 1 f f~L i i_ T ty Vi i | 

iTrKwrlfTi iHVTlrtNV ♦ NotK f trtTM NV ♦ crr~LN 


Ip-fttfB-XIbWAX l 

913 0FLTAV=0fcLTAVA2v© — 


0471th 



©1 -©477© 
~©t -04700- 





i»u fu ^iu ui usmiu 

914 CALL STINT ISOCt )tWtt)ATt:w^f Vttt 1 1 NOKI Kt , NHtnlNT ♦ NftmAXTI tM 0»81fr- 

rn 

uuul err — 

n^u 

Ir 

(Vtt— VBOAA) 9lb>904»413 01 XT4 

B30 

91 D Ufc 

t t AV*UtL 1 AV / 2 • 0 Or 04 

o«*tJ 


€0 TO 910 
904 VA* VAb 


XIPKLO-PKI 
GO TO 121 


200 HR 1Tb ( Jb*7009)0fcPTH 
KATIO * TOT lot/TUTOUT 


©1 04*60- 
HH- 04* 7© 


— Oi- 04090 
©1 04900- 


— Ol' 04920 
01 04930 



P0= (fc lN+tUUT ♦ /TU 
HKITt(JB*7011)fcIiX 
WKlTtl Jbt7012)fcUUT 
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602 entt STf-NT" (Tn>. 0 ,a. 0 ,mi)-, lcNGTU^t,! 
If- tPKLtJl 63O,630,6M 


-ortwofl- 
-01 oomo- 


— c all st int - ts 

Vt)=VtiO*HJOG1r 


■ »nbhaxt y~ 


01 0 5 2 40 


631 CALL MINT <SUC, -1.0, bflT t«H, VBU, 1 ,NGKi I 
- VtJ«VtiO * l-UUGt 


GO tt) 604 

603 X IL Im=X I BnAX — — 




MMISMB 


052YO 


~OT 00200 
Olr 06300 


621 X I L I rt=0. 6 

622 PlU=VU*(XlA-XILI )*Fltf-b 

CALL STINT ( SUC , 0. 0, bAl tmP , VbO, 1 , nGR I Rfc , Ntth I N I , ix H m A X T I 


01 05320 
01 05330 
01 05340 








Vb»VbO ~+~ FUOGfc 

CALL STINT <T,0.0,O.0*PKtO* 1 


[KJ 1 ^ <lt 




GO TU 604 

617 CALL STINT < SUC *-l .O* bATtrNf , VbO, l,NGKlPfc»N«NlN 

UhWH+ 4- mMtlHif 

TfNBHAXT) 


X hb»6PT0/ Vb I-< PKLO/T Vb-VOlGOfc I H 
GG TU 604 


- Vb-VbL * fOOGfc 
-*Tb«<PTG/Vbl-U 


0 544 0 



606 VU*VU+0. I 

—CALL ahps C 


GtF^TO 60^- 
60b CALL STIN T I SLHr r 


Oi 0553 0 

<\ | / vC C i i ' 

ul WV^XJ 


XldN- IbTO/VbT-fbKLO/T-Vb-VtHOOt H 
1F-H X IbN/X 1 b I— Or99T 619,619,618 


— Oi 0 5560 
— (Hr 05576 


GG TG 613 
619 Xlb-Xlb— 0.1 


6« IF<VdHAX-V6l 606.606*640 

620 CALL STINT ( SUL * X lb » bATfcNP, Vtt , 1 »NGK l t*fc,N8N INT-*NHmAXT 1 


- 01 0 56 00 

ft 1 i - C /■ 1 ft 


-01 05630- 


Xlb-IPTG/Vbl ~ IPKLO/I Vb-VOlOOtn— 

If TXILIh-XIbl 606,613*613 


-61 05650 
--Ol- 05660- 


I-H VO.LT * V€4 GO -TO 644 

O*** -p*HbT CUhPOT AT I ON CUMPUrTfc f UR- *>Nt‘ T I Mfc \ NCKfcHfcNT 


0-1 05670- 


PA»( VA+AOlUUfc )*XItt 
XlPKLO«PKLu/ < Vd-VOIUOt ) 


OT~- 05690— 

Ol 05700 


641 FLAG* 1.0 
GG TU 606 


700 TV5K=1000.0 

Vttb=AVnP5A-A0lUUt 
X I A = X I AM 


Ol 05730 
01 05740 
01 05750 


I 


70 





or 0576 a 


on 711 LKM= 1 , IMP ANfct 


01 05790 


711 CttNTlNUfc 


ai o5«oo 


CALL STINT (SUC,U.O,BATtNF,VBU, l,NGKlFt,NBNlNt 
VB=VBU*FUUGfc 


01 05 K 2 U 
01 05 H 30 


CALL STINT (T, 0 . 0 , 0 .U, 
vTu*ve 


^)1 05 K 50 
01 Tftntjrr 


1 NFWM, N INV » NCNV »NSfcK,tM-I NV * fc F 6 CNV , TN ) 
XlB=FTU/VB-FKLU/( VB-VUIUUt)-XlLT 


Ol 05HM0 


01 06 H 90 


701 CALL STINT (SUC,- 1 . 0 ,BAltMF,VBU,l,NGKIFfc,NBMlN 
VB= VBU*FUL)Gfc 


01 05910 
Ul 05920 


01 06970 
01 O 6 VH 0 


X 1 FKLU=FKLU/ ( VB-VUlUOt ) 
GU TlJ 704 


Ol 06000 
01 06010 


VTU-VB 

XIFKLO=FKLU/< VB-VOlUUt ) 


01 06030 
Ol 06040 


703 XlB»XlBWAX 

CaLL STINT ( SOC « X I B, BA I tnF, VB * 1 , rtGK I Ft * rNBnlttl , mBhAXT ) 


Ol 06060 
01 06070 


VT U* Vb+VKCRK 

CALL OK A I in (X IL,XILT ,FLO, VTU,CLSK, FVSK,tKSK,T ,inGK IFfc,iNSLT 


70 H XIn=FTO/VTU-XILf-FKLO/I VTU-VOlUUt 1 
16 <XlB-XlBmAX) 707 , 706,703 


Ol 06090 
Ol 06100 


CALL ANFS (XI A, VAB,F lfcWFS,X 5 ,XC,tlA,ANGLA,ANGLl fWGKlFC,NAUTU, 
lNCfcLLT ,AUlUUt»OfcLsr;G*KS»mUOUfcL*VUCCI , A DDT »UtLI I > 


01 06)00 
01 06190 
01 06200 


707 1 1 6 K = 1 0 

RCKK= ( VKCKK-OZCKK I /X I BNAX 
VTU=VBU+OZCKK 


1 NFWM,NINV,NCNV, NSfcK, tH- INV, tFFCNV, IN) 
X IB=FTU/VB-FKL 0 / ( VB-VUlUOt ) -X I LT 

_ 

01 

01 

05940 

03950 



GU TO 708 

706 IH VBMAX-VB ) 705 , 703,740 

.. . 

- -- - 01 06130 




















CALL OKAIn <XitfX-lLT*PLOf*TO*CLS*»TV5K*bKSPr1^N<HH"P*»i*5t-H 

1 NPWh* N InV ♦ NCnV ♦ NSbK * bPP I f*V * bPPCNV * 7 n I 

X lb -P 7U2 V TU - XH. T- P HLU2 4 V 7 It V OlUOU 

710 CALL STINT (SUC^XIb^bATbnP^Vtt^NGKIPbtNonlNTirNbinAXPT 
Vb-Vb*PUUGb 

VTU »VH « UgCHH * X i i n*6 fc KH 

CALL OK A I n 4XIL*XILT*PL0*VTu*CLSK»TYSK*bRSK *T **GK iPb »f*SL3-* 

INPWrt* N l NV * NCNV* NSbK , bPP | NV * bbtCNV * Tn I —— 

X *UN'-f TU/ V Tt l X 1L - T -P KLU/4 VI U -V D l UU Pl 

IF- (X.- -XI BN-O.Ol > 706 * 706 * 70 * 

70 * X I 0 * X i ON 

IT bH— I T bK- 1 

IP I 1 7 tK 1 7 Oft * 7 06 » 7 iO 
7*0 VC»V0*1.U 

I7-4~V A b .LT » VC > GO TtP 7 * 1 

70 * PA*XlA*< VAb+AUlOOt) 

X-1PKLO-PXLO/ < V fti-VOIOOt > 

V*-¥*rO 

60 TO 121 


61-0**10 

—01- 0**20 

64 - 0 *230 

- A 1 i ii, 3 i.i> 

Ol 0*260 

01 0 * 2* 0 

01 6*270 
Ol 6*2hO 

ol 002 * 0 

Ol 0*600 
Ol 0*610 

01 0 * 620 

- Ol 06660 


01 0*3*0 

Ol 0*3*0 

01 0 *3* 0 
Ol 06370 


7*1 


PLaG-I.O 
GO TO 7 O ft 


*00 OfttTb <6**101 NPAObL 

*10 PORMAl < IX * ‘NOnbbK OP PANbLS • 


• *13* »bXCbbOS NOorttK -A LLO Wb O 1 1 


6-P6P 

C**** ST 6 CUnPUT A f I On COnPLb*b POK UNb 
C SAI bLLITb bNbKGY bALANCb POKnA * S 
— 7000 P UH HA) — 1 31 2 ) 


T I Hb iNCKbnbNT 


01 

01 

o 4 - 


0*300 

063*0 

0 6 *00 


-7001 FOKrtAT <72H 


I 


7002 F UH HA T < 3 *H Uw Ab l te 70 HfcA U 4 A« L 6 * r ■ AOUH 7 - -i ■****- 


01 

Ol 

-64- 


06*10 

06*20 

0 6* 30 


7003 POKOAf -< 72 H 
700 * P OH ftA l < 12 73 4 


I 

-HOH- 


. 2 * 1 0-1 2 


01 

64 - 

1 7 7 73 0 4 


06**0 

06*60 

0 6*6 0 


N O . 1 3 * 1 H 7 H li t* f H I S O A < b U P 1 2 * 4 'H- 

H. J.FYLANO.H.KASHOSSbN POl 06*70 


-TOO* 


-1- SATbLLn* bLbCTPlCAL bNbKGY bALANCb 4C4 
2K0GK AnnbK S / 2 OHO lAttLb CUnntN 7 S « I 
P OK HA T — ( 1 8 H U - HU N CU nnbn3-6«> 


ill A 


7006 FUKn A I ( 30 HU NbW UK CHANGbO PAH Ant 7 bKS/ 1 H 0 ) 

7007 PUK«A 7 (I 3 .P 12 . 0 ) 

700 H P OK n A I (3 6H U P K H OK 4* T AB L E LOOKU P* N G H4PP — 1*4- 


(H 

Ol 

-64- 


06*00 

06610 

0 6* 20 


-P 1 6 r 3¥ P 7 


700 * FOKMAT ( 2 bH 0 UbPTH UP 01 SCHrtKGb 4 PbKCbN'l I P 7 . 3 / / > - 

7010 POKnA 7 < 10 H C /0 KAIIUP 6 . 3 //I 

■ Toll PO K N A 1 — ( * 3 H CH A HG b- bNb KGY ( 0 * 7 - 7 - nl N OT PS > 

7012 PUKi-iAT (* 3 H OISCHAKGb bNbKGY < HAT T-n 1 NU'fbS I -F 10 . 3//7 

7013 FOKnAT ( * 3 H UKH AVG PwK UISSIPATbU IN BAT TbR I bS 4 WATTS I F 10 . 3 / / I 

BUUl FUKHAl ( 1 7 **UH TnCNIGHT TIrtb) F 12 .*I 


01 06630 


01 06560 


1 ‘.C 


72 



r r - ■ «n .«,»>«»■ ' VWM - -iW<liil^» **»^ "••it* 9 •*' 



8002 

8003 


FORMAT 

FORMAT 


( I 7»40H 
( 17,400 


TU(0RBIT 
OELTA T., 


TIME) 


>12.4) 
.M2. A) 


01 06560 
01 06570 


80 04 FO RM AT -Ft 7 >40 0 V PM A X F 9 0LTS ) MM 1 


01 069 8 0 

01 06590 
01 06600 


8005 

8006 


FORMAT 

FORMAT 


( I7.40H 
< 17, 4 OH 


18MAX 

OZCRR 


(AMPS). 

(VOLTS), 


>12.4) 

>12.4) 


BOUT FO RMAT ttT ,4 PM VKC K* ( VU C TSt T; M l ? 1 . A ) 


0 1 06*10 

01 06620 
01 06630 


8008 

8U09 


FORMAT 

FORMAT 


< 17, 4UH 
(17 »4UH 


C/0 RATIO...., 
TVSR tVULTS), 


>12.4) 

> 12 . 4 ) 


8010 F O RM AT ( I T ♦ 4 0M t*SR ( W i n ! ) M l 2 .4 ) 


O f 0 664 0 

01 06650 
01 TI6660 


8011 

8012 


FORMAT 

FORMAT 


( I7,40H 
( I7,40H 


PIN. 

PLt). 


>12.4) 

>12>t 


8013 F OR M AT 1 t ' 7 ~ t 4 OM V U100 E ( VULT S ) Flg- ;4 t 


— 0 1 U 66 70 
Ol 06680 
— 01 06 690 


8014 

8015 


FORMAT 

FORMAT 


1 1 7 * 40H 
( 17,400 


BAMMAX tA-M). 
ETA-(OEG).. .. 


>12.4)- 
>12.4) 


8 016 F O RM AT ( 17 i40 H NS tT. ' . 1«t 


or-tto ro o 

01 06710 
Ol 06720 


8017 

8018 


FORMAT 

■FORMAT 


( 17,400 
ii 7 ,-40M 


18) 

+8> 


- TO 1 9 F O RM AT ( I7»40 O 


8020 

8021- 


FORMAT 

FORMAT 


07,400 
< I71-40H 


NBM inT i ... ,-rr. 

mbmaxT-. . . ...... w.. rwrv f.irnr.T; 

nCCLLT 

NPKt O. I nr 


>8t- 


- (H - -06 73 0 
<0t -867*0 

m .*■. t 1 » v 

trt On rot r 


-8 0 - 2 2 — F OR M AT (I 7»400 N t NV 10 ) 


01 0 6 7 6 0 


8023- 

8024 


FORMAT 

FORMAT 


( 17,400 
M7 t 


-Mtfty . . . . . . 


40H- -MSirR-.-.— 


-HH- 


-01—00770 
-OF -tJ 6700 


002 5 f OR M AT ( 1 7 ; 4 00 S VS TE m *t .Y ( ■ 1 .O - S mF T , O.O -n m , \ .tfM MM T )F } 2 . 6 ) 


-0 1 0 6 7 0 0 


8026 

8027 


FORMAT 

FORMAT 


f F7f40H 
(17-»40H 


NPRlNTv; 


FH-F 


■Ol — O6 80t ) 


HH- 


r\ l r ii l ji 

trt 1 J nn i tt 


0028 f OR M AT — ( ) 7 ,4 00 L ( ' » ' I ft V . 4. 12. 4 ) 


8020 

8031F 


FORMAT 

FORMAT 


t 17,400 
117,400 


ETFCMV. 

NBTEmP. 


-12-.AF 


1M 


Ol 0 68 20 
Oi 06 x30- 
-OF tt6840 


8031 F OR M AT (I7,40 H WOGRAU ; . . .1 8) 


8032 
803 3 


FORMAT 

FORMAT 


(17, 4 OH 
( 17,400 


8 03* F O RMA T — > t 7 t; 0 


8035 

8036 


FORMAT 

FORMAT 


( 17f40H 
(17, 400 


SIGISC (A/OtG Cl. 
SIGVOC (V/Ut« C). 

-of If PER C EN T) 

O I 208RCEMT )..... 
01 3 ( PERCENT ) 


...F12V5t- 
...F12.57 — 


0 1 0 68 00 
-01 06860 
01 06870 


, . r I g- . - 4 ) — 

, Ar) - 

;rM2-.4i— 


Ol 06 8 80 

-or 068V0 
01 06000 


8 037 FO RM AT ( 17,60 0 L) 1 6 ( MtrK C EM T ) F 1 2 . 60 


8038 

8039 


FORMAT 

FORMAT 


( 17,400 
( 17,600 


OVKPERCENT), 
0V2( PERCENT > 


>!?» 

>12.4) 


Ol 06 9 10 
-ot 06920 
01 06930 


8 040 FO RM AT ♦ 17, 6 88 1 w VPM O( V OLTS ) r r. 


8061 

8062 


FORMAT 

FORMAT 


{ 17,600 
( 17,600 


a i rmo ( amps ) . , 
AVOCO( VOLTS), 


, . . I 12.5) 

rr.-rre . -5-) 
,-r.F 12.tr) — 


0 -1 0 69 6 8 
01 06950 
01 069 60 


8 0 6 3 - F O RM AT 1 1 7 » 6 00 f Ot :T A ( P E R C EN T ) M2. 6 ) 


O l 06970- 
Ol 06980 
01 1)6990 
01 07000 


8066 

8065 

8066 


FORMAT 

FORMAT 

FORMAT 


( 17, 60H 
( I 7 , 6UH 
( 17, 60H 


TNUTdJEG C) 

no. OF CELLS/HATT, 
AUI0UE( VULT ) 


>12.6)- 

. 18 ) 

> 12 . 6 ) 
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-IttlT H 


8047 FORMAT +I7,40H PTtrH-vy.-. . . . . . .... 

8048 FORMAT U7,40M OtLTT (OtG €V 


8050 FuRMAT (I7.40M OtNFACtt 
8060 FORMAT f78H0 SOLAR ARRAY 1>rA 


1 O I L 'x I o i O • U i i ur Lr\« i 1 um t 

8V97 FORMAT < ?5M ^tASH-tft O NU T~ 




NS R=NGR I Ft 


TFR0F=O.-0 — 

-Hr=OvO 


A€-NV=0.0 

SER » 0 -.0~ 



XIt*PWML 

-1 IF TNlNV-il O f 3 t 4 


- 02-00230- 
-Of- 00240 


— Ft I f+StR — I I — 7 » 7 * 8 

CALL STINT 3T,O.O»0.O» AiNV»l»NGRII?fc»NlNV,NlNVI 


02 00270- 


OO TO 3 

6 CALL STINT IT tO.OtO.O f CNV, 1 » nGR I Pfc, nCNVtMOnV » 


Of- 00200— 
02 00300 


GO TO 5- 

8 CALL STINT (T tO»OyO*Ot StK » 1 1 NGR I Ft t NSER »NStR ) 


02 00320 
02 00300 


X I L = X I PROF - 
OIFF=VA— TVSR 
IF(L>IFF)34,34,35 


Of 00350 
02 00360 
02 00370 
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34 XILT=XIL+CLSR 
RET URN 


02 003H0 
02 00390 


36 XILT=XIL+CLSR+UIFF/ERSR 
RETURN 


02 00410 
0? 00420 


SUBRUUTlNEPRINTtNPRlNT,NPAGE,NTALLY,A,B*C»U,EtF*G,H,P,rj,R,StT,TJ»V)03 00010 
J A = b 03 00020 


I F ( NT ALLY >1*1*2 
l IF(NPAGE-bU)4,b,b 


03 00040 
03 000b 0 


WRITE ( JB » 9 ) 
NPAGt=0 


WRITE (JB»7)A»BfC»UtE»FtG»V*HfPt0»StT » R 
NTALLY=nPRInT-1 — - 


03 0007 0 
03 000H0 


03 00100 
03 00110 


2 N1'ALLY=NTALLY-1 

return 


03 00130 
0-3 00140 


6 FORMAT ( 120H1T Inb ACCOrr STAtfc- VTu VOInIGHTI tA- 03—001*0 


1 PA 


ISO -- I P K ttJ T EMP' 0 3 0017 O 


V FORMAT ( 43H 
END 


ST b UNLT>- VAwlDAYl/; 


lO IF 1Z-360.0I 2,2tl 

— 1- Z=Z-360.0 


03- 001 VO 

03 00200 


— 04-00020- 
04-00030- 


2 C«tL STINT < Z»u.O,U.U, a, 1,nGWIPB, 4,4>- 
ANGLt=A 


-04 00030- 

- -04 000*0 


tNt>— — — 04— OOOBO 

SOB-ROUT I Nt -ST AbH < SiGIS-ttSlGVUC»oISC,OV, fNOT.-KtrrtNOKi-F-ctNe-t-tLT-r- ~ltb OOOto 


01 MENS I UN VVEC ( 101 ) , XI VtC< lUl ) , TEMP ( 1 3 It VMAPl lOrlOlFt XIMAPl tbylOl +00 00031; 
DIMENSION VUCT ( lb ) tXI SCT ( lb) » VMP ( lbl» X IMPT 13 ) Ob 00040 


COMMON AttMP, AVMP, AXlMP,KEYt 
C NE* SUBROUTINE STASPt-TO UEGKAOt 


— — Ob OOObO 

G -FO“ NEW— PROCEEDOREOb 000*0 


JA=b 

Jb=6 

NGR I PE=U 


Ob OOOBO 

Ob 00090 
Ob 00100 
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' if (KtYti i,DOt«r,t 

I I F I KE Y ) 2000 1 3000,1000 


***************************** 

*jhiM 4 ufaiHiliiinto*iii ** mm *****<imi** 


€ IN FIRST cNtRY 



S1GVUC IS THtf A 
DISC IS AVRmO 


I- TUQT IS TNU I 

i KEY IS -1 


r. NCELLT IS NCCttT 

C— VUCOI IS VUCOE 


C UEtHF IS DELTT 

- 2 000 W E C UE-OfO- 


00-2001: l=2,iOi 

x evec < 1 1 =o. 1 1 


— oo zw m — ttttra 

CAtt. stint ivvEett»,o.i 


00230 

00240- 



2002 COnTInOE- 
GO TO 2005 


-+2W3 RtMN — 

-e EXTRAPOLATE INPUT CURVE 


- 2005 DU 2107 I - 2, lOt 




2107 CONTINUE 

- £106 SLDPE-I X I VEC ( w-1 I-X I VEC < N-2 > > / ( VVEC ( N-l 1-VVECIN-2 


PZCO = VVECHI - XlVtCU I/SLOPE - 

D« 2108 J ■ N» lOl - 


00 CURRENT DEGRADATION WITH GAMMA, TNUT , THETA 
X I SC"X I VEC I H 


DELT A I = < 1.0-GAmmA )*X I SC 
RECIP*1. O/GAMMA 
THET a = s igvuc 


05-0 02 0 0 - 

■ " /\ i\ *1 i > / ‘v 

XrD UXjduXj 


iV I f. f\ *, », 

05 TJtJjtU 


os oosoo - 


os-oosho- 

OS 00300- 


— OS 00410- 
-OS 00420- 


- OS — 00440— 
OS 00450— 


- OS 00470 - 
OS 00480- 


05 005 OO 
05 00510 - 


05 00530 
05 00540 
05 00550 






DELTAV=tTNt)T + 273. I6)*O.861At-0*F*ATTJG " CRtCt"P1 

05 

00560 

UELTAV ~ DfXTAV + ( 1.0 - T HE 7 A 1 *VUCt) I 

05 

00570 

DO 2600 I — 1,101 

09 

00980 

XIVtCMt* X1V2CII) - DttTAI 

05 

00590 

IP IXIVtet I ? t~ 2703, 2703,-2600 

05 

00600 


LOU, 101 


05 00620 
05 00630 


GO T 0-2006 

2600 CONT iNOt 


vw w v/wn I i iiw u 

; LOCATE POINT Ur 26*0 C 

— SLOPte = < X IVfcCtKtXt 


— 05 00650 
05 ^00660 


— 05 00680 
— 05- 00690 


otLTAv - (prco “ prcor-^-ottTAA 
-00-2 Tv* — I~2»lOl 


a 1 pmo=ov - 

AVPmU=0I St 


“00? to 

oo??o- 


I* WHIIIJ f»ll| 


05- 00? AO 


io "5 trCOrttJ C At t. 


NCE1_Lf-l5 r»etLLf 
V0C0I - IS OtNf-AC 


- -09 00900 


C OfcLTT IS OfcLT'f 
30UU ALl-A-SIGISC 


— 05 00920 
-09 00930 


FK I =01 SC 
7 K v = 0V 


05 00950 
09 00960 


UU 3001 I = 1,101 ----- - 

3001 VVfcC. ( I ) = VV tt ( I ) -X I VfcC ( i >*K£S 

StRItS K tS I ST AWC t OtGKAUAl I UN IS NUW LUMPLtFt 


05 00980 
05 U0990 
09 01000 








OeitKttil'tb VUC tOtftTb 1-+K5T 


00 3007 I » 2riOi 


iHX I VtC < I >13006,3006, 3007 


3007 CONTlNUt 


FlU. lbhPbRATORb VbCTOK 
IF <UbLTT) 3012,3012,3013 


— 30-1 3 CONTI NO b 


-VO€T I IT* VGC-AOOf*TKV 


-3202 V6CTI I )=VUC-<TbMP< I )-?NU> )* TKV 
: FKOANO I-V CORVt iNltT FAmI tY *r 


-PittbT TNNbK Otl 


ObtHFA V«T K V*HF 2 T 


X*rIT-TKl*7 2Tl*SiGVtK: - 
0 bNON * VOC— ObNF A<~T*V* 1 2 T 1 


03 OF3TO 


OMZ6- 1.0-26 


IF! OMZ6 >3300, 3203,3203 


— -3203 XfMAPT F,Jf*XIVfc£<GT-H(Kn*Gn£0- 
3040 -CONT iNOb 


00 4003 1-1,13 


03 013H0 


4000 PWKT=VMAPtl ,tO*XlHAP(I ,N) 


03 01410 


-1FIPWKT-PWR > 4002,4004,4001 


03 01430 


PWR=PWRT 
00 TO 4000 

4004 VMP< I )=VMAP( I ,IM) 


05 01440 
05 01450 


-6 

bbCUNO INNbH 00 F 


>3 Ole rO- 


Ot> 3203 0 — t,iOi 
Z-VMAtM If J>/UtNUM 


JJ UITnu 
lb. .iii ^ua 
























r 


■ 


X I MP ( I ) =X I MAP ( I »N1 

GU TO *005 “ ' 

*00? -S Ll t P E - lx i n A RI 1 » m - 1 ) -x I n a ? ( I ♦ * ) ) t (*p m 8 I | » n )- v w A Rt t ♦ w-tl-t- 
*003 VMP( I )=VMAP( 1 ,N-1 )+0. 00013 

XlMPI I )=XImAP< I,N~l)-StUPE*< VflPI I >-VmAP-M,i v-ltt 

Pw kT =v wp i 1 )»x 1 nr* < I ) 

IF (PWRT-PWR) *005, *005, *008 .... - 

*008 PWR=PWRT - — “ 

GU T fE * - 0 03 

*005 CONTINUE - . . 

*006 00 *007 1*1,15 — 

- -*0 07 - X fSCT I I ) - X I n AR (1 , 1) 

C EDIT RESULTS 

WRITE (00,7000) 

W RIT E ( JB,700l ) ( T EM8( J ) ,J-1,19 ) 

WRITE ( JB,7002XXlMP(J),J=l,15> - - — 

WRITE tt8,7O03t(VmR(tt, J=l,i5) 

W ft - ITL ( JB»700 6) m i CT ( I ) »l-l4l 9) 

WRITE- TtB,7007> ( XiSt?mrl=T, 19) 

00 “2 0*0 t=l-» 30 » 2 

W RITE ( J8 t 700* » JHXI -H f rRI 1 , J ) , 1 = 1 rl ~9) 

- 20 * 0 - hr tT E 1 10,7005 hi t tvintrtM ! , J t , 1=1 * 151 

WR I TE- T 10,7000 ) 

W RITL ( 18 1 700 I ) ( TL WPI J) i J-l , 15 ) 

WR i-TE i-10 ,7002 H- XIw RI It » 1 =7, 19 1 

WRtTE— <-10 r 70 0 3 1 ( Ami* H 1 > i=t-,t») 

W R l Tt f J 8 ,7006H9UCT ( I ) , I-t f 1 5) 

WRITE 1 18 *7007 1< XI SO I L I I t 1 =t ,~191 — 

OU 2 0*1-1 = 31,60,2 

W RITL (JO i 700 * ) J f < X Iw B R ( I , J > i t- 1 , 1 9 1 

-20*1 -WRITE T-tth7O05Tl, < VnAP< 1,1), 1=1 ,191 

WRITE 1 7000) 

W RITL <10,7001 ) ( TE MR ( J ) , 1 = 1 , 1 5) 

WR-ftE 118,7002 ) TX fwR < J ) ,J=1, 19) 

-WRITE ilB*7O03MVMPl 1)-,1=1*15> 

W RITL ( JB+7UU6) t VuC t t 1 ) » 1=1)1 9) 

WRITE < JB » 7007 )(XISCT(I)*!=1,19> “ 

00 20*2 J = 61,90,2 

W RITE ( J W »700* t l; <- X - lwA Ptl » J ) » I-l, 18 ) 

20*2 TtR ITE" ( J8, 7005 ) J » ( Vl*lAP( I , J) ,1 = 1, 15) 

WRITE < JB * 7000 ) 

W RITE ( J B , 700 1 ) I TLm R ( J ) » J — 1 » 15 ) 

WRITE ( Jb,7002> ( XIMPI 1) , J=l, 18) 

WRITE ( JB,7003)(VMP(1)»J=1,15) 

WRITE ( JB,7006) ( VUCTI I ) , 1=1, 15) 


0^-01*60 

05 01*70 

09 01* 8 0 

05 01*90 
— -09 01500 

05 01 - 9 t t) 

- 05 01520 

- 05 01530 

09 015*0 

05 01550 

-05 -0156U. 

0 5 01570 

03 01580 

05 11990 

0 9 0 t 6 0t ) 

05 01610 

05 01620 

09 01630 

05 01 6*0 

-85-01650 

0 5 01660 

05- Ot 670 

05 — 01 680- 

05 016 9 0 

05 0170 0 

0 5 01710 - 

0 5 01720 

05-01-7-30- 

05 017*8 

0 5 017 9 0 

85 0 17 6 0 

05 — 05770 

0 5 0)7 + 8 9 

05 01790 

05-01800 

05 01810 

05-01820 

05 01830 

0 5- 0 18*0- 

05 01850 

05 01860 

05 01870 

05 1)1880 

05 01890 

09 01900 
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WRITE -50tt»700?7l XI5CTI I 5»5*l» : 
00 2043 J * 51»84 f 2 


2043 WKITfc I JB,7005)J» I VMAPI I , J » » I « 1 ♦ 1 1> T 
RETURN 



rv/uv runnRi uirni cuir ur 

1 CfcLL I-V CUKVbS H)K 

icnrun« umc ucunwi 

tNtKbY rrALANLfc risrl 



ARG5 0L57O 
05 01 580 


— 701.12 FORMAT I14H0 I-MAX PWR 15F7.4) 

— 7003 FORMAT I 14H V-MAX PWR 15F7.4) 


05 02010 


7005 FORMAT I4H V ( 1 2 ♦ 8H ) 
-7006 FORMAT -( 14H VOLTS OC 


15F5.47— 
15F7.4) 



— C LOOK OF CURRENT toil) MAX FwR PtifiXT -G1VE15 VCtLt ^- 

iOOO - VCELi--=OISC 


/\ c iiOA/ n 

UzJ ViCVotr 


:MP<1) -HELMUT) /0trL5T 


NrH-NK+t- 


r~\ IT /\ r\ /\ 

tT5 uZutU 1 
- 05 - 02: 00- 



■ — -585)— 5334 » 533 3, 1354 

585 5NTO 5IG5SC FOT OPEN CIKCOH 




QrTl |Dai 

Kctunrl 

- 1334 IF I KEY — 588) 5 33 5 t 1 336 t 1338 


05 02 150 


-1530 005402 J - 2f l Ol 
N-»-i5 


05 021 8 0 - 
05 02550- 



- OO 1422 J - 2*101 
H • J 


05 02240- 
05 02250- 


1422 CONTINUE 

1423 CONTINUE 


1 VMAF I NH I ♦ N— 1 II2IXI MAP ( NH I » N ) — XI MAPI NH I f N-l > > 
GO TO 1210 


05 02270- 

- -05 02280- 

A t A A «»/W\ 


- 05 02300 

05 02350 


IFIVCELL-VMAPINLO, 1 ) >1100,1100, 1001 

1100 XILO-XlMAPlNLOt 1 ) 

GO TO 1010 


-05 02330 
05 02340 
05 02350 
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— riA=o.o or-oooso 

— OO lOtT LRM- 1 »NFANfcL 


IF (XCtLKMM ioo,iuu,i 07 OOOFO 

1 ANG=COS <fcTA/37.2Vb779:.)*CUS I ANGLA<LKm>/37r2*ttT7^>Si*r — 1 07 OOOVO 

-i nA/‘j7.2 * L » 77 93)» SI i' < A i nOLaI LK n > / 3 7. 2 <* 3 7 7 <»; > > » C« 5 ( t ANOtl t O? DO 100 

2LKl*t ) t-UtLANG) /57.29377V3) ** OT~OOtlO 

IF (ANG> 100» 100» 2 07 00120 



07 00140 
07 00150 
07 00160 



■A\r 



IT UOn-O.U- 
18 CUNT InUc 


1 Ai>OU/> 

UT UU3“v 

07-00390 


ACS * A 1)8 00040 

RETURN ON OOOSO 


4 A-O.O 

ACS - A 


IF (X+1.0I3»3 f 7 
A=»3. 14139263 


08 00070 
03 00080 


08 OOIOO 
08 001 10 


RETURN 08 00130 

6 A«aTan ( SORT ( 1.0-X*XmX> > 08 00140 


RETURN 

7 A«3»14139263+aTAN 
ACS = A 


(SORT (1.0-x*xmx»» 


08 00160 
08 00170 
08 00180 


100 CONTINUE 07— 00 4 i0 

1U1 RETURN 0T0042O 


FUNCTION ACS (X) 08 00010 

IF (X) 1,2»3 08 00020 


CALL STASH (UCV, VKNtttO.OfPTEHFS(LKN) »-tTi 
lAOOT »UfcLTT ) 


1VUC I ♦ AUU f »UELT T > 

UELV » UCV - 7 A I - OEL V 


THIRD ENTRY AN O LATtR ENT«ItS» TABLE L OOKU P 

CALL STaSH(TLU> VKNltO.OtPTENFSI LRN> »O.O t l«NORTPE t N C -E tt 


3010 DEL I a TLU*( 1 . 0-EFFEC f *CUG I 
PHIL*1.0/(tFFECI*CUG> 
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UPTURN 

OB 

00190 

fcND 

08 

00200 


U I MfcNS I UNNUMPT S ( 3 1 ) till 30) tt. 2(30) .13130) ,STG(2000 ) , DUMMY (30 >, AXft 20 6 00020 
1) — 0 6 00030 


fcOU I VALfcNCt (NAT. L3( 1 ) ) 
Ja = 3 


Oft 00090 
06 00060 


NGKlPfc=U 
I e < Kt Y ) 1.1.70 

-+-* 

NG= 1 

NUKmaL* 1 


06 00060 
Oft 60090 

-Oft-'O fttttft- 

06 001 10 

Oft 00120 


1397 PURrnAl(34Hl fnoLt 
GO TO 95 


CUNlfcTVfj T~ 


Oft 00140 
Oft 00190 


NOKhAL»2 
PRINT 1237 


06 OOl *6 
-Oft 604*6 


3000 KtTUKm 
779 WGK I Pt = l 


JO 66206- 

06-002 1 0- 


776 NGKlPfc=2 

WKITb < JB.9000) AKG1, ARG2,AKG3 ,I»i!n1bL .nAXTnt 


06 tnr 23 o 

wncw 


VUUO PUKWAT ( 20HO fcRKUK IN I LU» AKOl -M iT.3 »-6H AKOZ^M?r9 1 6 H 
IBM PllNT0t=14*8M MAXrtJl=14t 


BfcGlNNlNG UP STINT 
39 NUnfBt_=l 


riBTtTft 00260 
06 00270 


06 00200 

06-00300 


NUMt»TS< lt»0 
96 'P< I >69. 103, 102“ 


0ft00320 
06 00330 


37 PUKWaT < 2 A A * I A » 2 1 2 1 1 3 A4, A2 , i 2 ) 
PKlNl ll37»6fUAltfNAnt 


06 00360 


59 L8 = L1 (NUfiTBL ) 
N1-IL8-1 )/9+l 
00 68 IS-ltNl 


Oft 004 1U 

Oft 00620 
06 00630 



00 TO 62 
61 L6=L8 


L6=L3*NaT 

L7-l3*t4 


06 00 5 00 

iii* i |At \ IX 

XJTJ tAJTI t; 


L5=L2 ( rtOMl BL7 
tn*L5*l.6 


— 06 00530 
— 06 00040 


63 tftfi 65, 106 t 103 
103 KtAO 1 JA ,66 1 IOIJMmY (t0,t'*l,10)»lo 


107 STG< L3 ) •OOMmY ( 1 1 


-5ftH3l-U0MnYl*l 
65 K-K+l 


66 L 6 «t»omAT 
L7-LN-M.6 


16 1JJ-L41 67,68fOO 
67 JJ-JJ+1 


06 00560 
06 00570 



06 0 0650 
06 00660 


- 06 00600 
06 0065 O 



105 Lfcfc»NUM6TS<NUMTBL) + (L8*l )*U.5+1) 
16 ( L tt— 2000 1 llOO, 1101, 1101 


H02 NOiiP’f S< NUmTbL+ 1 ♦•tfcfc 

lOB NOtiT BL *nOMT bL + 1 


06 00740 
-06 00750 


06 00770 
r . v uV ho 


liol WKITfc < JB,llll)Ltfc 
GU TU 773 


06 00800 
06 00H10 


GU TU 775 

1111 6UKma1 < 17H TUU MANY 6U1WTS 181 


65 GU TU <773,776, 776), NG 

70 16 ( MlNTbL-MAXTBL 1 71,100*65 

71 UU 73 N AT=M InT BL , MAXT BL 


06 00830 
06 00840 


06 008 60 
06 00870 
06 00880 
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1-STG(K81) 
GU TO 9b 
91 IF (AhAXl 


>( ISEO>,STG< JT75TG(TSEO-l ) ,STG( J-l ) 1-9999. t9) 92 t 


06 01310 
06 Cl 320 
06 01330 
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92 FCT-< <STfi< U*ATfc>-ARG2m <STGt JJ >-AKGll*ST6< J-1T— 4 SF9FK91-A0919 06 01360 


2ARG1 )*STG ( IS tO) H / 4 4 STG4 I OAT£ 1-STG 4 A8 ) ) *4 bTG4 JJT-STG4*91 H 
GO TO 96 - 


06 01370 
06 013H0 


94 FCT=STG( !SfcO>-< STG4 JJ1-AKG1 >*4 STGH SfcO>-STG4 ISbO-1 1 T-/4ST64 J39- 06 01400 

1STG4K9H ■ - -- - 00-01410 


96 OOf*n¥4 H-FCT- 

0 7 Ml\T ■«9flT — 1 

7 r iTM I "T1H I X 


01960 


06-01660 










N I MB US -TfttttrtrS 
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